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Abstract 


A  finite  element  Investigation  was  conducted  to  analyze 
an  axial  tension  specimen  with  col linear  defects  placed 
symmetrically  about  a  center  crack.  The  material  modeled 
was  IN-718,  a  nickel -based  superalloy  used  in  Jet  engines. 
The  effects  of  crackxdefect  interaction  were  compared  using 
elastic,  elastic  plastic,  and  viscoplastic  constitutive 
models.  A  2-D  nonlinear  finite  element  code  called  SNAP  was 
used.  This  program  has  the  capability  to  simulate  crack 
growth  and  closure  by  releasing  or  closing  nodes  along  the 
crack  plane. 

Elastic  stress  intensity  solutions  were  developed  for 
two  different  finite  width  specimens.  The  stress  intensity 
versus  crack  length  plots  compared  well  with  infinite 
theory.  Results  reflect  the  defect  can  partially  shield  the 
crack  from  finite  width  effects.  A  critical  spacing  was 
also  noted  where  the  stress  intensity  of  the  crack  exceeded 
the  stress  intensity  for  the  combined  length  of  the  crack 
and  defect. 

Finite  element  analysis  of  a  crack/defect  configuration, 
considering  elastic-piastic  and  elastic -viscoplastic 
effects,  provided  crack  opening  profiles,  plastic  zone 
profiles,  and  stresS/^train  fields.  In  general,  the  defect 
has  a  prominent  influence  range  equal  to  approximately  one 
defect  length  for  all  constitutive  models.  The  presence  of 
a  defect  increases  the  magnitude  of  the  crack  opening  and 
stress>^straln  fields  in  front  of  the  crack  tip. 


I.  Introduction 


Background 

The  United  States  Air  Force  has  adopted  fracture 
mechanics  as  a  design  criteria  for  aircraft  engines  through 
the  Engine  Structural  Integrity  Program  CENSIP5  which  was 
instituted  as  a  Military  Standard  CID  in  1984.  The  damage 
tolerance  policy  is  based  on  "retirement  for  cause".  This 
policy  requires  the  determination  of  crack  growth  rates  such 
that  inspection  intervals  will  be  set  to  one-half  the  time 
required  for  an  existing  crack  to  grow  to  a  critical  size. 
Linear  elastic  fracture  mechanics  CLEFM5  provides  a  good 
estimate  of  remaining  life  for  components  subjected  to 
moderate  temperatures  and  loads.  However »  components  such 
as  turbine  disk  blades  are  subjected  to  both  high  stresses 
C200  ksi!)  and  high  temperatures  Cl 200  FD  which  leads  to 
interaction  between  creep  and  crack  growth.  For  these 
components,  the  Air  Force  Office  of  Scientific  Research 
CAFOSRD  has  set  forth  a  philosophy  that  damage  tolerances 
should  be  characterized  by  near-tip  stress,  strain,  and 
displacement  fields  with  emphasis  on  non-continuum 
description  of  structural  materials.  Application  of  damage 
tolerant  concepts  yielded  a  40  to  1  benefit  to  cost  ratio 
and  extended  part  life  by  a  factor  of  four  for  the 
FlOO-PW-200  core  engine  C25. 

Approach 

In  response  to  the  above  directives,  several  studies 
were  conducted  using  elastic -pi as tic  and  viscoplastic  finite 
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element  analysis  C3-9D.  The  geometries  considered  include 
both  axial  and  compact  tension  specimens  with  emphasis  on 
single  crack  phenomena.  The  goal  of  the  present  work  is  to 
build  upon  these  efforts  by  Incorporating  col linear  defects 
with  a  center -cracked  axial  tension  specimen. 

Two  primary  objectives  were  set  forth  in  this  study.  The 
first  objective  was  to  investigate  the  influence  of 
col linear  defects  on  the  stress  intensity  solutions  of  the 
center -cracked  tension  specimen.  This  was  accomplished 
using  the  J -Integral  method  developed  by  Rice  CIOD.  The 
second  objective  was  to  investigate  crack  growth  in  close 
proximity  to  the  defects  under  both  monotonic  and  cyclic 
loading.  The  material  was  1N-718>  a  nickel -based  suparalloy 
used  in  turbine  blades  for  Jet  engines.  The  material 
behavior  was  modeled  with  elastic-plastic  constitutive 
equations  using  the  Von  Mises  criterion  C113  and  the 
viscoplastic  flow  law  developed  by  Bodner  and  Partom  C1S> 
for  high  temperature  effects. 

Finite  element  analysis  was  conducted  with  a  program 
called  SNAP  C133.  The  SNAP  program  was  modified  by  Mercer 
to  allow  crack  growth  and  clraure  by  the  incorporation 
of  spring  type  boundary  nodes  along  the  crack  line  in  a 
manner  similar  to  Newman  CIO.  Additional  modifications 
were  made  to  allow  up  to  5000  degrees  of  freedom  and 
variable  crack  gro%d.h  rates  under  monotonic  loading. 
Puur-noded  linear  isoparametric  elements  were  used. 
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Literature  Review 


Single  Cracks 

As  previously  stated »  several  studies  have  been 
undertaken  using  elastic-plastic  and  viscoplastic  effects  to 
characterize  the  behavior  of  a  single  crack.  A  summary  of 
related  work  is  provided  here. 

Zahoor  and  Abou-Sayed  C3D  performed  elastic-plastic 
analysis  on  a  center -cracked  tension  specimen.  Both 
constant  strain  triangles  and  4-noded  quadrilateral  elements 
were  used.  Crack  tip  blunting  effects  were  noted  by 
prominent  element  rotation  at  the  crack  tip. 

Hinnerichs  CAD  examined  constant  load  creep  crack  growth 
for  a  center -cracked  specimen  of  IN-100  at  high  temperatures 
with  a  program  he  developed  called  "VISCO".  A  procedure, 
referred  to  as  the  hybrid  method,  was  developed  for 
determining  crack  extension  using  calculations  of 
viscoplastic  deformation  with  no  crack  growth.  In  this 
procedure,  the  difference  between  total  crack  deformation 
and  viscoplastic  deforoation  is  attributed  to  crack 
extension.  Extremely  good  crack  growth  predictions  were 
made  for  the  axial  tension  geometry. 

Nicholas  and  others  C6D  investigated  plasticity  induced 
closure  involving  short  cracks.  They  found,  from  an 
analytical  point  of  view,  that  closure  requires  some  amount 
of  crack  extension  in  order  to  develop  residual  strains 
behind  the  crack  tip.  The  plastic  wake  effect  developed 
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very  rapidly  for  fully  reversed  loading  but  required  longer 
propagation  distances  to  develop  under  positive  load  ratios. 

Wilson  and  Palazotto  C63  investigated  viscoplastic 
fatigue  in  an  IN-100  compact  tension  specimen  with  R  =  0.1. 
They  found  that  a  large  majority  of  plastic  straining  occurs 
within  the  first  three  load  cycles  and  the  stress  field 
remains  relatively  constant  after  one  to  three  cycles. 

Henkel  and  Palazotto  C7D  compared  viscoplastic  fatigue 
in  a  compact  tension  and  a  center -cracked  specimen  with  R  = 
-1.  They  found  the  size  of  the  plastically  strained  region 
at  the  crack  tip  is  a  major  factor  determining  the  amount  of 
closure  behind  the  crack  tip.  In  addition,  incomplete 
closure  behind  the  crack  tip  was  noted  at  full  negative  load 
for  both  specimens. 

Mercer  CSD  did  an  extensive  viscoplastic  study  of  a 
crack  growing  from  a  notch  under  cyclic  loading.  He  found 
that  the  notch  has  a  region  of  influence  equal  to  one  notch 
radius.  It  was  also  noted  that  a  highly  loaded  short  crack 
has  more  crack  tip  plasticity  than  a  long  crack  at  a  lower 
load  level  to  produce  the  same  stress  intensity.  Chestnut 
CO)  did  follow-on  research  for  large  round  notches  and 
proved  a  similar  influence  range  of  one  notch  radius. 

Multiple  Cracks 

Limited  publications  vwre  found  in  the  area  of  multiple 
crack  interaction.  The  work  summarized  below  discusses 
infinite  plate  analysis. 
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Matake  and  Imai  Cl 53  investigated  the  behavior  of  a 
small  col linear  defect  located  in  front  of  a  long  crack. 

The  possibility  of  pop-in  behavior  induced  by  the  main 
and  subcrack  was  demonstrated.  Pop-in  behavior  is 
characterized  as  an  abrupt  load  dropping  when  the  crack  and 
defect  combine  during  fracture  testing.  Under  constant  load 
this  phenomenon  is  characterized  as  an  abrupt  drop  in  the 
stress  Intensity  of  the  combined  crack  and  defect. 

Analytical  equations  for  the  stress  intensity  solution  at 
each  crack  tip  were  developed.  These  equations  are 
discussed  under  Fracture  Mechanics  Theory. 

Rose  C163  represented  microcracks  in  front  of  the  main 
crack  by  using  point -source  complex  potentials.  Simultaneous 
equations  for  parameters  characterizing  the  strength  of  the 
equivalent  point  sources  are  then  solved.  This  method  proved 
to  be  accurate  to  within  five  percent  of  the  exact  solution. 

Chang  C173  addressed  the  problem  of  noncoplanar  crack 
interaction  and  developed  equations  for  the  stress  intensity 
based  on  asymptotic  approximations  of  the  normal  stresses  at 
the  crack  tips. 

Yoda  Cie3  did  an  experimental  investigation  of  crack 
coalescence  in  glass.  He  noted  that  the  crack  velocity 
increases  with  the  stress  intensity  factor  as  the  cracks 
approached  each  other. 

Ang  C103  considered  two  eollinear  cracks  of  equal  length 
and  provides  a  boundary  integral  solution  to  this  problem. 
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He  also  presents  an  exact  solution  for  the  stress  intensity 
factor  for  the  outer  tip  of  the  cracks. 

Rubinstein  C203  developed  an  exact  solution  for  a 
macrocrack  interacting  with  a  periodically  distributed 
collinear  array  of  microcracks.  He  found  that  if  the 
microcracks  were  spaced  far  enough  apart  one  need  only 
consider  the  leading  microcrack  in  the  interaction  with  the 
main  crack.  For  this  case»  he  found  the  critical  spacing  to 
be  two  microcrack  lengths. 

Summary  of  Findings 

The  review  of  the  literature  reflects  significant 
accomplishments  in  2-D  analysis  of  single  cracks  under 
various  geometries  and  material  behavior  models.  Xn  the 
area  of  multiple  cracks,  substantial  work  has  been  done  for 
linear  elastic  analysis  in  an  Infinite  plane.  However, 
multiple  cracks  under  finite  geometries  and  nonlinear 
effects  is  completely  open  for  study. 
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II .  Theory 


Equations  of  Motion 

Finite  element  modeling  is  a  well  developed  technique 
for  approximating  the  “exact**  equations  of  motion  for  a 
body  undergoing  deformation.  The  SNAP  finite  element  code 
developed  by  Brockman  CISD  is  based  on  an  incremental 
virtual  work  expression  that  characterizes  a  body  deforming 
from  a  state  Ck~ll>  to  a  state  CkD  and  is  given  by: 


=  I  dVol  *  J  -  J 


Ck-11  Ca.  15 


Vfhere: 


^ijkl  Iricremental  materiel  constitutive  matrix 
e|^^  B  Incremental  linear  strain 


'ij 


■  Incremental  nonlinear  strain 


B  Surface  Tractions 


*  Body  Forces 

e  2*'**  Piola-ICirkhoff  stress  tensor 
6C  ^  B  Virtual  quantity 

Mercer  C83  provides  a  detailed  derivation  of  the  above 
equation. 

The  corresponding  finite  element  expression  in  total 
Lagrangian  form  is  written  Cneglecting  body  forces^  as: 


ki«t>  *  “=9>)  <“>  •  k<’»  -  k-l«> 


C2.23 


7 


Wher«:  =  Tangent  stiffness  matrix 

CK  3  =  Geometric  stiffness  matrix 

g 

<u>  -  Nodal  displacement  vector 

<T>  =  External  force  vector 

<I>  s  Internal  force  vector 
The  total  stiffness  matrix  is  divided  into  tangent  and 
geometric  components  to  account  for  nonlinear  material 
behavior  and  large  displacements  respectively.  The 
nonlinear  plasticity  problems  were  solved  using  the  initial 
stiffness  method.  This  method  uses  the  stiffness  matrix 
formulated  with  the  first  increment  of  displacement  for  the 
entire  solution.  During  a  load  increment  the  internal 
forces  in  the  elements  may  not  be  in  exact  equilibrium  with 
the  applied  forces.  This  extra  force  is  considered  a 
residual  force  and  is  applied  to  the  next  time  increment. 
Equilibrium  for  the  current  load  step  is  met  when  the 
residual  force  and  displacement  corrections  are  less  than 
user  specified  tolerances.  A  detailed  description  of  the 
initial  stiffness  technique  is  given  by  Owen  and  Hinton 
Cll>. 

Constitutive  Models 

Three  constitutive  models  were  used  in  this  analysis  for 
calculating  stresses  and  strains  in  the  elements.  A  brief 
description  of  each  is  provided  below. 


d 


Linear  Elastic 


The  linear  elastic  constitutive  model  follows  the 
classical  relation: 

<o>  =  tD'i<cy  ca.  35 

T 

where:  <o>  -  ^  ®’l2^ 

egg  ej^g> 

[Dl  =  Elastic  constitutive  matrix 
El astic-Pl astic 

For  this  model,  total  strains  are  decompx)sed  into 
elastic  and  plastic  parts: 


e 


ij 


ca.45 


Incremental  plastic  strains  are  calculated  via  a  Prandtl* 
Reuss  relation  C115  as: 

de^P  *  dX.  Ca.B5 

Where:  dX  •  3  do 

a  e  H 


®1J  "  D»''t»«-orle  .trM.  - 

'  ^  " 

H  «  Slope  of  stress-plastic  strain  curve 
o  »  Effective  stress  «  Jj  1 3/2CS.  .S.  ,53 


The  corresponding  matrix  form  of  the  equation  is  given  by 
Yamada  Cai3: 

<do>  «  C  DepKde>  C2. 65 

where:  do  »  increawNnial  stress 

de  e  incremental  strain 


O 


CDep}  =  ttlastic-plastic  constitutive  Matrix 
Mercer  C83  shows  the  elastic -plastic  matrix  in  expanded 
form. 

For  this  model,  a  Von  Mises  yield  criterion  is  used  to 
determine  when  the  material  exceeds  a  critical  value  of 
recoverable  elastic  energy.  For  the  case  of  uniaxial 
tension,  yielding  occurs  when  the  effective  stress  exceeds 
the  yield  stress  of  the  material. 


Bodner-Partom  Viscoplastieity 

The  Bodner-Partom  flow  law  accounts  for  viscoplastic 
behavior  as  well  as  rate  sensitivity  and  strain  hardening 
effects  C1B3.  The  model  description  given  by  Mercer  C8D  is 
presented  here  for  completeness. 

For  small  strains,  the  total  strain  rate  is  decomposed 
into  elastic  and  plastic  parts  by  : 


*lj  “  ""ij  ■**  ®1J 


C2.75 


Vfttere:  ^ij  Elastic  strain  rate 

«  Viscoplastic  strain  rate 

The  viseoplastic  strain  rates  are  calculated  using  a 
Prandtl -Ruess  type  relation: 


where  is  the  deviatoric  stress  and  is  the  second 
deviatoric  stress  invariant. 

The  constant  n  controls  the  model's  strain  rate 
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sensitivity,  while  Dq  is  the  maxi mum  value  of  strain  rate  in 
shear.  The  Z  parameter  is  an  internal  state  variable  which 
accounts  for  the  degree  of  material  work  hardening  and  is 
expressed  as: 


2  =  Z^  +  CZq-  Z^^exp  C-m  Ca.  OD 

where  Z^  and  2^  are  the  material's  initial  and  maximum 

values  of  hardness  respectively,  and  the  constant  CnD 

controls  the  rate  of  work  hardening.  The  term  W  accounts 

P 

for  the  plastic  work  including  thermal  recovery  of  hardening 
at  high  temperature  and  is  defined  as: 

9p  =  J ♦  J  CS.10J 


The  rate  of  thermal  recovery  of  hardness  Is: 

.r 


rec 


■^{^1 


ca,ii> 


where  Zg  is  the  minimum  e)cpected  value  of  the  hardening  at  a 
given  temperature.  The  constants  CiO  and  CrD  are  material 
parameters  «^ich  are  chosen  to  match  low  strain  rate 
^secondary  creep)  test  data.  Secondary  creep  is  defined  as 
the  balanced  condition  when  the  rate  of  work  hardening 
eguals  the  rate  of  thermal  recovery  C4).  the  recovery  term 
becomes  essential  In  high  teo^Merature  analysis. 

Stresses  are  calculated  using  the  following  equation: 

<do>  «  toi  t<ds>  -  cds^i  ca.ia) 


%duNre  tDl  is  the  elastic  stiffness  matrix. 

Owen  afid  Hinton  Cll)  set  forth  the  following  equation  for 
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C3.133 


calculating  tho  plastic  strain  rate: 

Where  for: 

^  =  0  Euler  integration  scheme  Cfully  explicits 

^  s  1  Fully  implicit  scheme 
1 

^  g  Crank -Ni  col  son  rule,  or  semi -implicit 
A  semi -implicit  scheme  is  en^loyed  in  the  SNAP  code  to 
calculate  the  viscoplastic  strains.  This  scheme  is 
unconditionally  stable  as  described  by  Hughes  and  Taylor 
CS2D  which  allov^  larger  load  increments  to  be  applied 
during  the  solution.  However,  accuracy  is  decreased  If  the 
load  steps  become  too  large.  Chestnut  CO)  provides  a 
discussion  on  how  the  viscoplastic  solution  proceeds  in  the 
SNAP  code. 

Beaman  C83)  determined  the  Bodner-Pariom  material 
parameters  for  XN-7ie  listed  in  table  2.1  below. 


Table  2.1:  Bodner-Partom  Parameters  for  XN-718  •  1200  F 


Parameter 

Description 

Value 

Elastic  Modulus 

23.6  X  10  ksi 

n 

Strain  rate  exponent 

3.0 

% 

Strain  rate  limit 

10®  s«c“‘ 

Initial  hardness 

23S.3  ksi 

Maximum  hardness 

260.3  ksi 

Minimum  hardness 

104.1  ksi 

m 

Hardening  rats  exponsnt 

2.875  ksi 

A 

Hardsning  rscovery  coeff 

1.6  X  lO'^sc”^ 

r 

Hardening  recovery  exponent 

7.0 
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Linear  Elastic  Fracture  Mechanics 


Fracture  mechanics  studies  the  effects  of  externally 

applied  loads  and  specimen  geometry  on  crack  tip  stresses* 

displacements*  and  growth  rates.  ^aecimen  thickness  and  the 

stress  intensity  factor  provide  the  similitude  p»arameters 

necessary  for  linear  fracture  toughness  comparisons.  Broek 

C243  states  that  the  analysis  is  considered  linear  as  long 

as  the  plastic  zone  in  front  of  the  crack  tip  is  small* 

which  is  when  the  stress  is  low  with  respect  to  the  yield 

stress  Co  <  0. 8  o  When  the  stresses  are  higher*  the 

ys 

plastic  zone  will  spread  beyond  the  point  at  which  it  is  a 
unique  function  of  the  stress  intensity  factor.  The 
development  of  fracture  aechanies  is  well  described  by  Broek 
C243  and  therefore  will  not  be  elaborated  upon  here.  The 
equations  used  for  this  analysis  are  briefly  described 
below. 

Center "Cracked  Plate 

For  linear  analysis*  the  stress  at  the  crack  tip  is 
expressed  by  the  stress  intensity  factor  which  for  atode  one 
loading  Copening  mode)  is  given  by: 

B  /?  o  yi*a“  C2.i4) 

Where:  ft  «>  finite  boundary  correction  factor 

o  «  far  field  stress 
a  »  half  crack  length 

Xn  order  to  apply  equation  2.14  to  a  finite  element 
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model »  the  proper  correction  factor  for  the  finite  width  and 
height  was  required.  Isida  C25^  developed  finite  boundary 
corrections  for  a  variety  of  boundary  conditions  with 

the  use  of  Laurent  expansions  of  con^lex  potentials 
satisfying  the  stress  free  relations  along  the  crack.  Table 
2.2  lists  a  representative  example  of  the  finite  width 
corrections  for  an  axial  tension  specimen  subjected  to 
uniaxial  loading.  The  accuracy  of  the  finite  boundary 
corrections  are  regarded  as  correct  up  to  four  figures  C2SD. 
The  half  height/half  width  ratios  of  0.4  and  1.0  were  used 
in  the  mesh  validation  efforts  of  this  study. 


Table  2.2  Finite  Boundary  Corrections  C2S:> 


cA> 

0.4 

0.6 

0.8 

1.0 

0.0 

1.0000 

1.0000 

1.0000 

1.0000 

0.1 

1.060 

1.033 

1.021 

1.014 

0.2 

1.266 

1.130 

1.083 

1.065 

0.3 

1.620 

1.286 

1.184 

1.123 

0.4 

1.843 

1 . 407 

1.323 

1.216 

0*6 

2.247 

1.773 

1.406 

1.334 

0.6 

2.806 

2.123 

1.702 

1.481 

0.7 

3.670 

2.650 

1.040 

1.680 

Where:  y  *  crack  half  lengtlv>^peci men  half  width 

*  specioMMn  half  height/specimen  half  width 

Macrocrack  -  Nultiple  Defect  Clnflnite  Flate> 

Ihe  presence  of  defects  or  large  inclusions  near  primary 
cracks  in  ductile  materials  have  been  observed  in 
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experimental  work.  It  is  known  that  the  growth  and 
nucleation  of  such  defects  in  the  vicinity  of  a  primary 
crack  play  an  important  role  in  high  temperature  creep 
growth  C4^.  The  presence  of  defects  can  also  shield  the 
primary  crack  thereby  directly  affecting  the  toughness  of 
the  material  C16I>.  However »  quantitative  effects  of 
crackxdefect  interaction  are  not  very  well  understood.  The 
aim  of  this  work  is  to  provide  some  insight  into  this 
interaction. 

For  this  analysis,  a  defect  is  defined  as  a  through-the-* 
thickness  crack  that  is  approximately  one-fifth  or  less  than 
the  length  of  the  primary  or  macrocrack.  The  substantially 
smaller  defect  length  makes  a  favorable  argument  for 
limiting  crack  growth  to  the  macrocrack  tip  only. 

A  crack  with  col 11 near  defects  placed  symmetrically 
about  the  y  axis  as  shown  in  Figure  2.  i  was  selected. 

Col linear  analysis  simplified  the  finite  element  mesh 
development  and  and  allowed  direct  comparison  to  existing 
analytical  equations  for  the  stress  intensity  factor  as 
described  below. 


-c 


c  b  a 


SL 


4- 


LCl  4  2S  4-  8S6^ 


♦ 


Fig.  2.1  Macrocrack  -  Multiple  Defect  Geometry 
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Matake  and  Imai  Cl 5!}  studied  the  case  of  a  tnacrccrack 
with  two  col linear  defects  for  an  infinite  plate.  They 
obtained  analytical  expressions  for  the  stress  intensities 
at  each  tip  of  the  crack  and  defect  for  the  case  of  a 
uniform  tensile  stress  normal  to  the  crack  plane.  The  mode 
one  stress  intensity  solutions  for  the  above  geometry  are: 


K  =  o 
a 


‘^b  =  - 


^  Cl  +  as  +  as63 
a 


nL  Cl  aSdi 

a 


1  /  1  .  ECIO  1 

k  1  ^  KCk3  I 

n  r  1  ECk3  1 
k  1  a  KCk3  ~  J 


ca.  15D 


ca.  163 


K  = 
c 


o  nL 

4  2 


nL  i  ECkD 

a  n  Rnr3 


Ca.l73 


....  ,  1  ♦  s  4  as6 


^S5* 


ca.i83 


ca.  103 


n  «  Jci  -  k®3 

o  B  far  field  stress 

5  B  defect  ratio  of  crack  length 

6  s  spacing  ratio  of  defect  length 

KCk3  and  ECk3  are  complete  elliptic  integrals  of  the 
first  and  second  kinds  and  are  given  by  the  following 
equations: 

d4> 


KCk3  a 


I 


n 

g 


ca.ao3 


•Ici  - 


k^sin^<K) 
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ECkD 


Ca.  313 


n 

Jci  -  k^sin®^3 

0 

Stress  Iniensiiy  Methods 

The  elastic  stress  intensity  factor  can  be  determined 
for  crack  problems  In  finite  element  analysis  using  the 
J-lntegral.  compliance  methods,  and  extrapolation  of  the 
stress  and  displacement  fields  near  the  crack  tip.  For  this 

I 

study,  the  J -Integral  and  stress  extrapolation  technique 
were  used.  A  description  of  each  method  Is  given  below. 

J  -  Integral 

Rice  Cl 03  has  shown  the  J -Integral  as  defined  along  a 
contour  around  the  crack  tip  CPig.  3.23  Is  the  change  in 
potential  energy  for  a  virtual  crack  extension  Cda3, 


Fig.  3.2  J-Xntegral  Contour  Path 
The  J-lntegral  is  given  by: 


Wh®r®; 


j  =  J  (  Wdy  -  ds] 


C2.22^ 


W 


=  strain  anergy  dansity 
=  CElasticD 


mn 


C  El  as ti  c  -PI  asti 


*  traction  vector  along  contour  path 
u  s  di spl acamant  vector  along  contour  path 
da  s  increment  of  distance  along  the  contour 
s  alastic  strain  tensor 

fPj  =  plastic  strain  tensor 
»  stress  tensor 


The  stress  intensity  factor  is  then  calculated  from  the 
J-integral  by: 

K^=  /T  J  CPlane  Stress?  C2. 23? 

Where:  E  »  Modulus  of  Elasticity. 

Stress  Extrapolation 

With  this  technique*  the  stresses  normal  to  the  crack 
tip  are  used  to  determine  the  stress  Intensity  factor.  The 
stress  in  the  vicinity  of  the  crack  front  is  given  by: 

a  m  K.  cosC^>^  Cl  +  sinC^>/2?  sinC^^l  C2.24? 

J' 

V"i2nr' 
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For  the  case  of  ^  =  O*  the  resulting  expression  for 

K,  is: 

1 

K,  =  y  2nr  a  C2.2SD 

1  y 

Where:  o  -  stress  in  y  direction 

r  =  radial  distance  from  crack  tip 
The  averaged  stress  intensity  value  for  the  crack  tip  is 
extrapolated  to  r  «  o  based  on  the  method  described  by  Broek 
C20. 


Irwin  Plastic  2Sone 

The  basic  theory  behind  linear  elastic  analysis  of  a 
crack  states  that  an  infinite  stress  field  is  present  at  the 
crack  tip.  In  actuality,  this  cannot  occur  since  the  metal 
will  plastically  deform  once  the  yield  stress  is  met.  The 
plastically  deformed  region  is  known  as  the  plastic  zone  and 
a  rough  estimate  of  its  size  was  set  forth  by  6.  Irwin  as 
discussed  by  Broek  C24>.  The  approximation  is  a  circular 
region  in  front  of  the  crack  tip  with  a  radius  given  as: 


C2.26) 


ys 


Where:  R_  »  plastic  zone  radius 

P 

a  •  far  field  stress 

naterial  yield  stress 
ys 

a  «  crack  half  length 

A  comparison  bet%NMn  the  above  equation  and  the  plastic  zone 
size  generated  with  the  finite  element  technique  follows  In 
the  discussion  of  results. 
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III.  Pinito  Elttmant  Models 


This  analysis  modeled  a  0. 1  inch  thick  axial  tension 
specimen  subjected  to  a  uniaxial  load  normal  to  the  crack 
plane.  With  this  approach,  the  effects  of  finite  boundary 
conditions  normally  present  in  experimental  work  will  be 
analytically  modeled.  Figure  3.1  depicts  the  configuration 
of  the  crack  and  defects  in  the  axial  tension  specimen.  Due 
to  symmetry  only  one  fourth  of  the  specimen  was  modeled  with 
the  finite  element  mesh. 

Gener at ion 

Two  different  size  meshes  were  generated  to  evaluate  the 
finite  width  effects  on  the  crack/defect  configuration.  The 
first  mesh  CPig.  3. is  a  symmetric  mesh  with  a  height  and 
width  of  0.5  inch.  The  second  mesh  CPig.  3.33  has  a  height 
of  0.  S  inch  and  an  extended  width  of  1.S3S  inches.  Pour 
noded  isoparametric  elements  were  placed  in  a  reduction 
pattern  to  transition  smoothly  into  the  small  elements  along 
the  crack/'defect  region.  The  crack/'defeci  region  of  both 
meshes  consists  of  0.002  inch  square  elements.  This  size 
was  selected  to  correspond  with  the  meshes  generated  by 
Mercer  C83.  Mercer  proved  this  size  of  element  in  the 
crack  tip  region  gave  adequate  stress/stral n  data,  plasticity 
solutions,  and  stress  intensity  solutions  for  modeling  crack 
growth  with  a  similar  magnitude  of  stress. 
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I 
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riGUHE  3.2  tVNNEYHlC  NOOIL  OF  2H12L  TDISION  tPEClNXN 


The  crack /'defect  region  for  this  analysis  was 
approximately  five  times  larger  than  the  work  of  Mercer  C83 
and  Chestnut  COD  and  the  number  of  elements  required  to 
generate  the  models  increased  on  the  same  order  of 
magnitude.  This  raised  a  concern  that  the  computer 
processing  time  for  nonlinear  analysis  would  become 
prohibitive.  To  reduce  this  problem,  the  aspect  ratio  of 
the  elements  was  increased  to  a  range  between  two  and  four 
as  the  mesh  diverged  from  the  crack/defect  region.  Cook 
C26D  states  that  aspect  ratios  up  to  two  will  yield  good 
stress  results  and  aspect  ratios  as  high  as  seven  will  yield 
good  displacement  results. 

The  program  PDA-PATRAN  was  used  to  develop  the  finite 
element  meshes  needed  for  this  analysis.  Due  to  the  size  of 
the  meshes  C2070  and  3288  D. O.  P.  respective! yD  it  was 
beneficial  to  band  the  resulting  stiffness  matrix  to 
conserve  memory  and  processing  time.  A  matrix  is  "banded" 
if  all  nonzero  coefficients  cluster  about  the  diagonal. 

Cook  C80D  provides  detailed  information  about  optimization 
techniques.  A  program  called  SWAP  CAppendix  8D  was  used  to 
convert  the  neutral  files  generated  by  PDA-PATRAN  to  NASIRAN 
input  decks  to  aecoa^lish  the  bandwidth  ^timization. 
Optimization  reduced  the  bandwidth  from  28d  to  76.  The  SWAP 
program  was  then  used  to  convert  the  optimized  neutral  file 
to  a  SNAP  input  deck. 
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Validation 


To  validato  tho  mosh  configuration  fo\jr  oxamplas  of  a 
center  cracked  plate  were  performed  on  the  symmetric  mesh. 
The  far  field  stress  was  12  ksi  normal  to  the  crack  plane. 
Both  the  J-integral  and  stress  extrapolation  technique  were 
used  to  determine  the  stress  intensity  for  each  crack  size. 
Figure  3.4  shows  an  exan^le  of  the  stress  extrapolation 
technique  using  equation  2.24.  The  averaged  stress 
intensity  value  for  the  crack  tip  is  extrapolated  to  r  ^  o 
based  on  the  method  described  by  Broek  C&4^.  Equation  2.14 
with  the  corresponding  finite  width  correction  given  by 
Isida  C25^  for  the  particular  crack  size  was  used  to 
calculate  the  th€K3retical  value  of  the  stress  intensity. 

The  results  were  accurate  to  within  8  percent  or  less. 
Tables  3.1  and  3.2  provide  a  summary  of  the  crack  sizes  and 
results. 

Table  3.1:  Center  Cracked  Plate  Stress  Intensity 


Crack 

a 

1C 

**1 

Theory 

K 

^1 

Stress 

Percent 

Errw 

0.100  in 

7.10 

6.06 

1.80 

0.144  in 

Q.03 

8.80 

1.50 

0.176  in 

10.30 

10.30 

1.20 

0.200  in 

11.57 

11.68 

1.00 
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Table  3.2:  Center  Cracked  Plate  Stress  Intensity 


I 

I 


Figure  3.5  depicts  the  path  location  used  for  the 
J-lntegral.  Three  different  length  paths  were  used  with 
negligible  differences  Cl ess  than  0.1  percents  being 
recorded  for  paths  ranging  from  30  to  04  elements.  The 
elements  in  the  paths  have  an  aspect  ratio  of  one  and  a  size 
of  0.002  in.  This  analysis  and  the  work  of  Mercer  C8d 
confirm  the  independence  of  the  J -integral  routine  to 
elenwnt  size.  However*  it  was  found  that  the  algorithm 
contained  in  the  SNAP  code  requires  the  element  aspect  ratio 
be  less  than  two.  As  long  as  the  path  selected  followed 
this  criterion*  the  J-integral  results  were  Independent  of 
path  length.  Feu  aspect  ratios  between  two  and  four*  up  to 
seven  percent  error  below  the  theoretical  values  was  noted 
in  the  calculations  for  the  stress  intensity.  The  writer 
believes  this  error  is  accumulated  in  the  weight  function 
technique  used  for  the  contour  path.  This  technique 
extrapolates  the  stress  values  at  the  eleoMit  gaussian 
points  C containing  the  contour  path3  to  points  along  the 
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FIGURE  3.4  STRESS  EXTRAPOLATION  OF  STRESS  INTENSITY  FACTOR 


MU 

flOIRE  4  *  IHItCllftl.  MM 
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path  itself.  Those  elements  with  aspect  ratios  greater  than 
two  slightly  underestimate  the  true  stress  field  as  noted  by 
Cook  Cc!6>.  The  small  error  in  each  element  is  then 
accumulated  along  the  contour  path  resulting  in  a  total 
error  of  seven  percent. 
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IV.  Results  and  Discussion 

The  problem  of  a  crack  interacting  with  col linear  defects 
was  analyzed  with  linear  elastic,  elastic-plastic,  and 
viscoplastic  constitutive  equations.  The  results  of  each 
constitutive  model  are  presented  below. 

Linear  Elastic  Analysis 

The  purpose  of  the  linear  analysis  was  to  compare  the 
the  finite  element  solutions  of  crack /defect  interaction  to 
existing  infinite  plate  theory.  The  results  given  below 
show  that  finite  element  analysis  provided  comparable 
solutions. 

The  first  objective  was  to  determine  the  influence  range 
of  a  defect  on  the  crack  stress  Intensity  factor.  Three 
different  length  defects  were  studied  with  the  symmetric 
mesh.  As  pointed  out  in  section  IX »  a  defect  is  defined  as 
a  through-the-thickness  crack  that  is  approximately 
one-fifth  or  less  the  length  of  the  primary  crack.  The 
J -integral  method  was  used  to  measure  the  stress  intensity 
of  a  crack  subjected  to  a  far  field  load  of  16  kips  €36  ksi^ 
normal  to  the  crack  plane.  Rubenstein  €303  used  the  method 
of  complex  stress  potentials  to  analyze  the  defect  influence 
range.  He  determined  the  defect  influenced  the  stress 
intensity  of  the  primary  crack  at  a  distance  equal  to  tvK» 
times  the  defect  length  €8a3.  Based  on  his  efforts*  the 
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Initial  crack /def«ct  spacing  was  set  to  two  times  the  defect 
length  for  this  study.  At  this  di stance »  the  addition  of  a 
defect  provided  a  one  percent  or  less  increase  in  the  stress 
intensity  for  a  center  cracked  axial  tension  specimen.  The 
crack  tip  was  then  moved  toward  the  defect  in  0.  002  inch 
increments  until  a  three  percent  Increase  was  present  in  the 
stress  intensity  over  the  center  crack  solution.  As  shown 
in  table  4.1.  a  distance  equal  to  1.1  defect  lengths  was 
required.  Once  inside  this  range,  a  sharp  increase  in  the 
stress  intensity  gradient  begins  to  occur  which  is  also 
reflected  in  the  analytical  solutions  in  Figure  4.1. 


Table  4.1:  Defect  Influence  on  Stress  Intensity  Solution 


Crack 
a  Cin3 

Defect 
a  Cin3 

Spacing 

Cin3 

CCP 

ksiVin 

CCP  w^Defect 
ksiVin 

Percent 

Increase 

0.214 

0.010 

0.022 

38.47 

37.47 

2.77 

0.104 

0.015 

0.032 

33.18 

34.32 

3.44 

0.172 

0.020 

0.  044 

30.67 

31.73 

3.70 

M  Far  Field 

Stress  ■ 

:  38  ksi 

An  analytical  comparison  was  then  conducted  to  determine 
if  the  defect  Influence  range  on  the  crack  stress  Intensity 
was  independent  of  crack  length.  Figure  4.1  shows  a 
comparison  between  the  center  cracked  plate  solution  Ceq. 
2.143  and  the  macrocraek -multi pie  defect  solution  Ceq.  2.173 
for  crack  sixes  raivging  from  0.132  -  2.016  inches,  the 
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DEFECT  SPACING  (A)  (INCHES) 

STRESS  INTENSITY  VERSUS  CRACK  SIZE  COMPARISON 


results  indicate  that  crack  length  doesn’t  have  a  measurable 
effect  on  the  defect  influence  range.  Furthermore*  the 
analytical  solution  confirmed  a  prominent  defect  influence 
range  of  approximately  one  defect  length. 

The  second  objective  in  the  linear  analysis  was  to 
analyze  finite  boundary  influences  on  the  crack  and  defect 
stress  intensity  solutions  and  compare  the  finite  element 
results  to  the  macrocrack -multi pie  defect  solutions  given  by 
equations  2.  IS  -  2.17. 

Figures  4.2  -  4.4  show  a  comparison  of  the  crack  stress 
intensity  versus  the  defect  sp»acing  for  the  same  three 
defects  used  in  the  first  objective.  A  far  field  load  of  6 
kips  C12  ksiD  was  applied  to  the  specimen.  The  6  kip  load 
ensured  the  size  of  the  plastic  9one  in  front  of  the  crack 
tip  would  be  insignificant  for  the  lincrar  analysis.  The 
stress  intensity  results  via  finite  element  analysis  compare 
<well  with  the  plots  given  by  infinite  theory.  Furthermore* 
a  small  four  to  five  percent  decrease  in  the  stress 
intensity  range  was  noted  for  the  extended  width  mesh 
Cfigure  3, 3D  which  is  1.5  times  wider  than  the  symmetric 
mesh  Cfigure  3. 2D.  The  small  drop  in  the  stress  intensity 
indicates  that  the  a  defect  will  essentially  shield  the 
crack  tip  from  the  finite  width  of  the  plate.  The  majority 
of  the  increase  in  the  crack  stress  intensity  over  infinite 
theory  appears  to  be  due  to  the  finite  height  of  each  mesh. 
The  finite  height  in  each  case  is  constraining  the  area  over 
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which  the  stress  field  can  act  thereby  forcing  a  direct 
stress  increase  in  the  local  region. 

An  individual  comparison  of  Figures  4.2  -  4.4  shows  that 
as  the  defect  increases  in  size,  so  too  will  the  influence 
range  of  the  defect  and  the  magnitudes  of  the  stress 
Intensity  solutions  for  a  given  crack/defect  spacing.  For 
this  analysis,  doubling  the  defect  size  resulted  in  a  ten 
percent  increase  in  the  stress  intensity  at  a  crack/defect 
spacing  of  2  mils. 

Figures  4.2  -  4.4  also  show  that  the  crack  stress 
intensity  will  exceed  the  stress  intensity  of  the  combined 
length  C crack  defects  when  the  defect  is  in  close 
proximity  to  the  crack.  This  occurs  since  two  stress 
singularities  Ci.e.  crack  and  defect  tips!)  in  close 
proximity  will  generate  a  larger  stress  concentration  at  the 
primary  crack. 

Figures  4.5  and  4.6  provide  a  direct  analysis  of  the 
spacing  where  the  crack  stress  intensity  will  exceed  the  the 
stress  intensity  of  the  combined  length  Ccrack  *  defect)  for 
the  symmetric  and  extended  width  mesh  respectively.  This 
spacing  is  termed  a  '*critical  spacing"  by  Tamake  and  Xmai 
Ci8).  If  the  crack  tip  is  inside  the  critical  spacing  and 
the  far  field  stress  is  increased  such  that  the  crack  stress 
intensity  equals  the  toughness  of  the  material  Ci.e.  stress 
Intensity  for  catastrophic  crack  growth)*  the  crack  will 
propagate  into  the  defect.  When  the  crack  combines  with  the 
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STRESS  INTENSITY  VERSUS  DEFECT  SPACING 


CRACK  -  DEFECT  SPACING  (A)  (INCHES  M  10  -3) 

FIGURE  4.3  CRACK  STRESS  INTENSITY  VERSUS  DEFECT  SPACING 


FIGURE  4.4  CRACK  STRESS  INTENSITY  VERSUS  DEFECT  SPACING 


CRACK  -  DEFECT  SPACING  (A)  (INCHES  M  10  >3) 

FIGURE  4.5  CRITICAL  SPACING  FOR  SYMMETRIC  MESH.  P  -  6  KIPS 
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FIGURE  4.6  CflXTICAL  SPACING  FOR  EXTENDED  HESH,  P  -  6  KIPS 


d«f<»ct,  tho  str«ss  intansily  will  drop  bolow  ih«  toughnoss 
value  to  the  the  combined  crack  and  defect  stress  intensity. 
An  increase  in  load  will  then  be  required  to  increase  the 
new  main  crack  stress  intensity  to  the  toughness  value  for 
propagation  to  continue.  This  phenomenon  is  referred  to  by 
Tamake  and  Imai  Cl 53  as  pop-in  behavior. 

The  critical  spacing  for  the  defects  used  in  this 
analysis  ranged  betvwen  6  to  8  mils  which  is  approximately 
40  to  60  percent  of  the  respective  defect  lengths.  Doubling 
the  defect  size  resulted  in  a  33  percent  increase  in  the 
critical  spacing  for  both  the  symmetric  and  extended  %d.dth 
meshes.  Here  we  note  that  defect  size  definitely  increases 
the  critical  spacing*  but  due  to  the  limited  data  no  exact 
correlations  bet%MMn  defect  size  and  critical  spacing  can  be 
made.  A  five  percent  drop  in  the  critical  spacing  %ias  also 
evident  between  the  symmetric  and  extended  width  mesh 
indicating  finite  width  has  a  small  influence  on  the 
critical  spacing.  The  small  influence  of  the  width  is  again 
due  to  the  defect  shielding  the  primary  crack. 

Figures  4.7  **  4.0  show  a  comparison  of  the  stress 
intensity  of  the  outside  tip  of  the  defect  to  infinite 
theory.  Again*  comparative  trends  to  infinite  theory  %«ere 
obtained.  For  this  case*  the  stress  intensities  calculated 
with  the  extended  width  mesh  converged  very  close  to 
infinite  theory.  This  result  follows  the  intuitive 
prediction  that  the  outer  tip  of  the  defect  will  be 


FIGURE  4.7  DEFECT  STfl^  INTENSITY  VERSUS  DEFECT  SPACING 
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FIGURE  4.0  OEFECT  STRESS  INTENSITY  VERSUS  DEFECT  SPACING 


CRACK  -  OBJECT  SPACiNB  (A)  (INCHES  M  iO  ~3) 
ricM  4.9  DEFECT  STRESS  INTENSITY  VERSUS  DEFECT  SPACING 


primarily  influ*nc«d  by  width  rathor  than  the  height  of 
the  specimen. 

The  linear  analysis  conducted  so  far  in  this  study 
validated  the  program  and  finite  element  meshes  against 
available  theory. 

Elastic -Plastic  Analysis 

The  elastic -plastic  analysis  used  small  displacement 
formulation  for  calculating  strain  displacement 
relationships.  The  range  of  the  problem  considered  was 
established  by  the  nonlinear  validation  efforts  performed  by 
Mercer  C8D  with  the  SNAP  code.  He  determined  that  small 
displacement  theory  was  valid  for  an  elastic  stress 
intensity  of  3S  -  45  ksiVin  for  a  center  cracked  plate  with 
a  material  yield  stress  of  140  ksi.  Based  on  his  efforts* 
this  study  selected  35  ksiVin  as  the  elastic  stress  ^ 

intensity  to  be  used  for  initial  crack  'lengths^  l1his..v 
al  1  owed  an  i nvesti  <^ti of  crack  growth  wi thi  n  a  ipd asti c  - 
regicw  as  'well  as  lmi^y  a  strain  smII  enough  so  that  small 
displacement  theory  could  be  Imposed. 

To  limit/the  amount  of  computer  ti«M*  two  primary  cracks 
of  initial  half  length  Ca>  B  0.106  inch  and  Ca:>  e  0.^05 
inch;  and  one  defect  of  half  length  Ca>  «*  O.Ol  inch  were  ' 
selected.  In  addition*  the  eraek/defect  spacing  was  set 
that  the  plastic  sones  ahead  of  the  crack ^and  defect 
overlapped.  The  propt^  load  Ci^  and  10  Klps^  was  then 
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appli«d  to  th«  spocimon  to  gonorato  a  stress  Intensity  of  35 
kslVln  for  each  crack  length.  At  full  load,  the  element  in 
front  of  the  crack  tip  had  a  stress  of  160  ksi  and  six 
percent  strain. 

This  portion  of  the  analysis  used  the  elastic -plastic 
relations  of  Eqs.  2. 4  -  2.6.  The  crack  could  not  be  grovm 
in  this  part  of  the  analysis  since  atten^ts  to  grow  the 
crack  toward  the  defect  resulted  in  negative  values  along 
the  diagonal  of  the  stiffness  matrix  during  equation 
solution.  This  was  occurring  since  the  stress<^strain  curve 
is  essentially  flat  at  six  percent  strains.  As  a  result, 
the  elastic-plastic  analysis  consisted  of  individual  cases 
of  crack  size  and  defect  spacing. 

Figures  4.10  and  4.11  show  the  shape  of  the  plastic  zone 
developed  for  a  fixed  defect  position  and  increasing  crack 
lengths.  Xn  Figure  4.10  the  primary  crack  is  approximately 
half  the  size  of  the  primary  crack  in  Figure  4.11  but  is 
subjected  to  a  higher  load  to  develop  the  same  stress 
intensity.  We  note  that  the  plastic  zones  for  the  two 
figures  are  approximately  equal  indicating  that  similitude 
with  the  crack  stress  intensity  gave  similitude  in  the  size 
of  plastic  zone.  As  the  free  surface  of  the  crack  was 
placed  closer  to  the  defect,  the  size  of  the  plastic  zone 
decreased.  An  analogy  to  this  effect  is  to  consider  the 
strain  energy  present  In  the  crackxdefect  spacing  as  water 
flowing  through  the  area  between  the  crack  and  defect.  As 
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the  area  becomes  smaller*  the  amount  of  water  flowing  or  In 
this  case  the  magnitude  strain  energy  present  in  the 
crack/defect  region  is  restricted. 

Equation  2. 26  is  compared  on  the  first  plots  of  Figures 
4.10  -  4.11.  This  calculation  confirms  the  crack  and  defect 
must  be  in  the  given  ranges  in  order  for  the  plastic  zones 
of  the  crack  and  defect  to  overlap. 

Figures  4.12  and  4.13  show  the  corresponding  crack 
opening  profiles.  It  is  seen  that  the  elastic-plastic 
analysis  generates  a  blunting  of  the  crack  tip  ^Ich  was 
essentially  constant  for  all  cases.  Based  on  the  pattern 
that  developed  with  the  given  profiles*  it  appears  that  one 
can  estimate  what  size  crack  will  not  have  bliinting 
characteristics  by  extrapolating  a  line  to  the  crack  plane 
as  shown  in  Figure  4.12. 

Figure  4.14  shows  a  direct  comparison  between  an  elastic 
and  an  elastic -plastic  crack  opening  profile.  This 
comparison  highlights  the  blunting  characteristic  in  the 
elastic-plastic  crack  opening.  Furthermore*  plasticity 
effects  decreased  the  strain  energy  present  at  the  crack  tip 
as  measured  by  the  J-integral  technique.  As  shown  in  Figure 
4.14  this  decrease  yielded  a  68  percent  drop  in  the  crack 
stress  intensity  factor  for  the  given  geometry  and  load 
condition. 
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LOCATION  (INCHES) 

FIGURE  4.14  CRACK  OPENINS  CONPAflISON.  P  -  25  KIPS.  100%  LOAD 


Figure  4.15  illust.ra'tes  the  influence  of  the  defect  on 
the  crack  opening  profile.  A  16  percent  increase  in  opening 
was  evident  along  the  entire  crack  profile  for  a 
crack/defect  spacing  of  0.026  inch. 

The  strain  profiles  betv^een  the  crack  and  defect  are 
shown  in  Figures  4.16  and  4.17.  For  both  primary  cracks* 
there  is  an  increasing  strain  profile  as  the  crack/defect 
spacing  is  decreased.  In  each  figure*  we  note  the  strain 
curves  are  concave  up  which  is  a  direct  result  of  increasing 
strains  at  the  crack  and  defect  tips. 

In  looking  at  the  magnitude  of  the  strains  present  for 
the  various  crack/^defect  spaclngs*  it  becomes  questionable 
as  to  how  much  strain  can  exist  in  the  plastic  zone  of  the 
material  before  instantaneous  failure  at  load  would  occur. 
Since  no  experimental  data  is  available  for  this  work*  one 
must  attempt  a  comparison  to  existing  finite  element  work  by 
others.  One  assumption  is  to  consider  failure  would  occur 
when  strains  along  the  crack  plane  approach  six  to  eight 
percent.  This  was  the  strain  level  Mercer  C8)  required  to 
match  experimental  crack  growth  to  his  model  of  a  compact 
tension  specimen.  This  is  a  rather  global  assumption  since 
we  are  making  a  comparison  between  an  axial  tension  specimen 
and  a  compact  tension  specimen.  The  only  similitude  between 
this  work  and  Mercer's  is  the  size  of  the  elements  along  the 
crack  plane*  the  material  thickness*  and  the  material  itself 
CIN-7183.  We  note  the  size  of  the  elements  used  in  a  finite 
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FI6UBE  4.i5  CRACK  0PENIN6  CONPARISON.  P  -  18  KIPS.  lOOX  LOAD 


CRACK  —  8  -  0.106  -  0.130  IN 
DEFECT  -  a  -  0.010  IN 
P  -  25  KIPS 
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DIStANCE  FfWM  CRACK  TIP  (INCHES  *  10  -3) 

F16URE  4.16  ELASriC/PLASTIC  STRAIN  PROFILE  BETWEEN  CRACK  AND  DEFECT 


OISTANCe  FROM  CRACK  TI»  (INCHES  «  10  -3) 

FIGURE  4.17  ELASTIC/PUSTIC  STRAIN  PROFILE  BETWEEN  CRACK  AND  DEFECT 


olement  modol  will  hav«  a  direct  effect  on  the  magnitude  of 
the  strains  and  therefore*  plays  an  inqaortant  role  vAien 
conqaaring  strains  between  models.  As  stated  earlier.  Mercer 
also  used  elements  sized  to  0.002  inches  along  the  crack 
region  and  since  no  other  work  was  available  with  this  size 
element;  his  work  was  selected  for  a  comparison.  If  one 
uses  this  criterion  based  on  the  above  clarifications*  we 
can  see  from  the  figures  that  failure  at  approximately  10 
mils  in  Figure  4.  IG  and  12  mils  in  Figure  4.17.  By 
selecting  a  maximum  strain  criteria  for  failure*  one  can 
postulate  that  plasticity  will  increase  the  critical  spacing 
given  by  the  stress  intensity  factor  in  the  linear  analysis. 

Viscoplastic  Analysis 

Viscoplastic  modeling  at  1200  was  accomplished  with 
the  Bodner~Par tom  algorithm  in  the  SNAP  program,  the  time* 
temperature*  and  rate  dependent  material  behavioi*  was 
accounted  fior  with  experimentally  determined  coefficients 
listed  in  Table  2.1. 

To  allow  comparison  with  the  elastic -aplastic  analysis* 
the  initial  crack  half  length  and  defect  half  length  were 
set  to  Ca3  «  0.2CG  inch  and  Ca!)  w  0.01  inch  respectively, 
the  same  external  load  or  IB  kips  C3G  ksi3  was  also  used  for 
the  viscoplastic  analysis. 

Crack  growth  under  both  monotonic  and  cyclic  load  was 
considered.  Crack  growth  is  accomplished  in  the  SNAP 


program  by  the  addition  of  "springs"  to  the  boundary  nodes 
along  the  crack  plane.  The  boundary  condition  of  each  node 
is  set  by  the  spring  stiffness.  Free  nodes  have  zero 
stiffness  and  constrained  nodes  have  a  stiffness  equal  to 
1000  times  the  maximum  diagonal  entry  in  the  structure 
stiffness  matrix.  User  defined  criterion  set  the  node 
release  tolerance  and  growth  rate  of  the  crack  as  a  function 
of  the  load  cycle.  A  detailed  description  of  this  method  is 
given  by  Mercer  C85. 

The  monotonic  and  cyclic  load  cases  are  depicted  in 
Figure  4.18.  For  the  monotonic  load  case»  a  fixed  crack 
growth  rate  Cda/'dtD  of  0.004  inch/^second  was  selected.  For 
the  cyclic  load  case,  a  fixed  crack  growth  rate  per  cycle 
Cda^dnD  of  0.004  inch,^cycle  was  selected.  However,  no 
restriction  was  placed  on  the  growth  rate  Cda/'dt)  during  the 
load  phases  of  the  cyclic  case.  The  above  rates  were 
selected  to  limit  the  amount  of  computer  processing  time. 

In  actuality,  thousands  of  cycles  are  required  for 
cyclic  crack  growth.  To  attempt  such  an  analysis  with 
finite  elements  would  be  economically  impossible.  However, 
Wilson  and  Palazotto  C6:>  state  that  most  of  the  plastic 
strain  properties  of  fatigue  loading  are  characterized  in 
one  to  three  cycles.  Thus,  the  limited  number  of  cycles 
used  in  finite  element  analysis  can  still  yield  a  good 
approximation  to  the  plasticity  induced  features  present  in 
the  problem. 
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Figure  4.18  Viscoplastic  Load  Ratio  Cases 


Figures  4.10  -  4. 8R  depict  the  viscoplastic  zone  shapes 
as  the  crack  is  grown  toward  the  defect  for  the  monotonic 
and  cyclic  load  cases.  For  the  initial  crack/^defect 
configuration,  the  height  of  plastic  zone  is  approximately 
20  mils  for  monotonic  loading  and  16  mils  for  cyclic  loading 
as  compared  to  6  mils  for  the  elastic*^lastic  analysis. 

Here  we  see  that  applicati«xi  of  the  Bodner-Partom 
viscoplastic  flow  law  at  high  temperature  in  combination 
with  crack  growth  has  increased  the:  size  ef  the  plastic  zone 
height  by  a  factcx*  of  2-'3  depending  on  the  type  of  load 
used. 

For  both  load  cases,  we  also  see  that  the  free  surface 
of  the  defect  constrains  the  initial  plastic  zone  from 
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riCURE  4.81  VlfCOPUISTIC  SQNl  SHAPES  f  CYCLIC  UMO 
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riCURE  4.28  VltCOPUISTlC  SONl  SHAPES  /  CVCtlC  UMD 


propagating  forward  with  the  crack  tip.  This  constraint 
caused  a  smaller  "super  strained"  region  to  be  formed  inside 
the  initial  strain  hardened  region  as  the  crack  moved  toward 
the  defect.  The  "super  strained"  region  is  the  smaller  area 
inside  the  original  CshadedD  strain  hardened  region  as 
depicted  in  the  figures.  Here,  "super  strained"  means  the 
remaining  region  that  still  has  stresses  exceeding  the  yield 
stress  and  a  corresp>onding  increase  in  the  strains  over  the 
previous  crack/defect  configuration.  As  noted  with  the 
elastic-plastic  analysis,  the  size  of  the  "super -strained" 
region  continues  to  decrease  as  the  crack  approaches  the 
defect . 

Figure  4.23  shows  the  plastic  zones  generated  between 
the  crack  and  defect  with  crack  growth  and  compares  these 
zones  to  a  fixed  crack  of  the  same  dimension.  The  increased 
profiles  reflect  the  time  and  load  dependent  nature  of 
plasticity.  A  15  to  50  percent  increase  in  the  plastic  zone 
area  was  due  directly  to  crack  growth  in  the  viscoplastic 
analysis. 

The  crack  opening  profiles  for  crack  growth  under 
monotonic  and  maximum  cyclic  load  are  shown  in  Figures  4.24 
and  4.25.  Crack  growth  under  monotonic  load  CFig.  4.24D 
resulted  in  a  sharper  crack  tips  as  compared  to  the  maximum 
cyclic  load  CFig.  4.253  which  shows  blunted  crack  tips. 

Again,  the  amount  of  blunting  was  essentially  constant  for 
all  the  crack/’defect  sp>acings  considered. 
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rXQiSS  4.33  UXSCOPUISTIC  ZONE  SHAPE  COMPARXSQN 


LOCATION  (INCHES) 

FIGURE  4.24  VISCOPUSTIC  CRACK  OPENING  PROFILES. 
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LOCATION  (INCHES) 

FI6URE  4.25  VISCOPLASTIC  CRACK  OPENING  PROFILES.  R  -  0.1.  100%  LOAD 


Both  load  cases  reflect  the  presence  of  a  plastic  wake 
forming  over  the  crack  as  it  moves  through  the  strain 
hardened  region.  This  presence  is  evidenced  by  the  abrupt 
drop  in  the  crack  opening  profile  near  the  crack  origin  and 
by  the  Jagged  contour  of  the  diminished  profiles.  The 
plastic  wake  formed  immediately  under  monotonic  crack  growth 
but  formed  one  cycle  CO. 004  inch^  later  for  the  cyclic  load 
case.  The  decrease  in  the  crack  opening  profile  due  to 
residual  plastic  deformation  is  highlighted  in  Figure  4. 26 
by  the  fixed  crack  of  equal  length. 

Figure  4.27  shows  the  crack  profiles  at  minimum  cyclic 
load.  Crack  closure  is  established  within  one  cycle  and  the 
amount  of  closure  doubles  as  the  crack  moves  toward  the 
defect.  Careful  attention  must  be  given  to  the  markers  used 
to  identify  the  individual  crack  profiles  to  see  the  closure 
Czero  displacements  for  a  particular  crack  length.  Figure 
4.28  highlights  the  closure  phenomenon  of  one  crack  profile 
by  comparing  it  to  an  elastic  solution. 

A  comparison  of  a  viscoplastic  crack  opening  profile  to 
an  elastic  profile  at  full  cyclic  load  is  presented  in 
Figure  4.20.  Again*  the  presence  of  the  plasticity  reduces 
the  strain  energy  at  the  crack  tip  with  an  associated  44 
percent  decrease  in  the  stress  intensity  factor.  The  stress 
intensity  was  calculated  with  the  J-'integral  and  application 
of  equation  2.23. 
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FIGURE  4.28  VISCORLASTIC  CRACK  OPENING  PROFILE  COMPARISON 


LOCATION  CINCHES) 

FISmE  4.27  VISCOPLASTIC  CRACK  OPENINS  PROFILES.  R  -  0.1.  lOX  LOAD 


LOCATION  (INCHES  «  10  -3) 

FISUflE  4.28  e^ACK  OPENINS  PROFILE  COMPARISON 


FIGURE  4.29  CRACK  OPENING  PROFILE  COMPARISON 


For  both  monotonic  CFig.4. 30D  and  maximum  cyclic  load 
CFig.4.31I>»  a  constant  stress  profile  was  reached  along  the 
crack  plane  when  the  crack  had  grown  within  one  defect 
length  Ci.e.  0. OS  inchD.  The  cyclic  analysis  stabilized  to 
a  slightly  higher  stress  of  180  ksi  versus  175  ksi  for  the 
monotonic  case.  For  lO  percent  cyclic  load  CFig  4.3S>»  a 
fully  compressive  field  is  generated  for  crack/'defect 
spacings  ranging  from  S  to  12  mils.  Between  IS  and  28  mils» 
we  see  a  stress  field  that  starts  compress! ve»  goes  positive 
and  then  decreases  as  it  approaches  the  defect.  When  the 
defect  is  20  mils  or  one  defect  length  from  the  crack  tip  we 
see  that  the  stress  at  the  defect  is  zero. 

The  strain  profiles  for  monotonic  and  cyclic  loading  are 
shown  in  Figures  4.33  -  4.35.  In  Figures  4.33  and  4.34.  we 
see  the  same  concave  up  strain  profiles  as  noted  in  the 
elastic-plastic  analysis.  However,  for  the  case  of  10 
percent  cyclic  load  in  Figure  4.35,  we  see  a  combined 
concave  up  and  concave  down  shape  in  the  strain  profiles. 

The  inflection  points  along  these  curves  correspond  to  the 
location  where  the  compressive  stress  field  converts  to  a 
positive  tensile  field. 

As  with  the  elastic-plastic  analysis,  the  question  of  a 
critical  strain  level  must  be  addressed.  If  one  uses  the 
failure  criteria  of  eight  percent  strain  for  the 
crack/defect  spacing,  the  possibility  of  instantaneous 
failure  occurs  at  approximately  10-12  mils  which  is  the  same 
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DISTANCE  FROM  CRACK  TIP  (INCHES  ft  10  -3) 

FIGURE  4.31  VISCOPLASTIC  STRESS  PROFILE  BETWEEN  CRACK  AND  DEFECT 


eoh>(cin'vmcM«4oo)GDrsujio^mcu^ 
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FIGURE  4.33  VISCOPLASTIC  STRAIN  PROFILE  BETWEEN  CRACK  AND  DEFECT 
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FieUf^  4.34  VISCOPLASTIC  STRAIN  PROFILE  BETWEEN  CRACK  AND  DEFECT 
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DISTANCE  FROM  CRACK  TIP  (INCHES  *  10  -3) 

FI6URE  4.35  VISCOPLASTIC  S^IN  PROFILE  BETWEEN  CRAC3<  AND  DEFECT 


spacing  given  in  the  elastic -plastic  analysis. 

A  second  and  far  more  realistic  failure  criterion  in  the 
cyclic  analysis  is  the  growth  rate  Cda/dtD  during  a  load 
increment  of  a  cycle.  As  stated  earlier,  the  growth  rate 
per  cycle  Cda/dnD  was  set  to  0.004  inch/cycle.  To  meet  this 
criterion,  two  nodes  were  released  at  a  user  specified 
stress  of  160  ksi  in  each  cycle  of  loading.  ^ecifying  the 
stress  for  node  release  allowed  the  node  release  time  to 
be  an  independent  variable  in  the  solution.  The 
corresponding  release  times  were  then  recorded  by  the  SNAP 
program.  By  subtracting  the  release  time  of  the  first  node 
from  the  release  time  of  the  second  node,  the  exact  time 
increment  for  a  0.004  inch  growth  is  known  and  the  growth 
rate  in  a  load  cycle  Cda.^dtD  can  be  calculated.  Table  4.2 
shows  the  growth  rate  results.  We  see  the  growth  rate  is 
initially  0.016  inch/’sec  and  then  decreases  to  a  constant 
rate  of  0.0125  inch/sec  for  the  next  four  cycles.  For  the 
seventh  and  eighth  cycles,  the  growth  rate  is  infinite.  By 
selecting  the  sixth  cycle  as  the  last  cycle  for  stable 
growth,  we  have  a  crack/'defect  spacing  of  0. 006  inch  lor 
instantaneous  fracture.  This  is  a  33  percent  increase  over 
the  0.006  inch  prediction  given  by  the  linear  stress  intensity 
predictions.  Using  this  method  for  predicting  failure 
yields  a  slightly  larger  strain  profile  of  approximately  ten 
percent  in  the  crack^^defect  spacing  as  shown  in  Figure  4.  34. 
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Table  4.2  Cyclic  Crack  Growth  Rate  Cda/dt:) 


Cycl  e 

dax’di 

Cycle 

da/dt 

Hf 

C i n/sec^ 

# 

C i n/sec5 

1 

we 

5 

0.0125 

2 

0.018 

6 

0.  0125 

3 

0.0138 

7 

Infinite 

4 

0.0126 

8 

Infinite 

Initial  Loadi^  No  Growth 

4n  additional  check  was  conducted  to  see  if  the  primary 
crack  length  would  affect  the  influence  range  of  the  defect 
when  working  with  viscoplastic  flow  law.  Figures  4.36  " 
4.38  show  the  stress-strain  profiles  a  crack  that  is 
approximately  half  the  length  of  the  original  viscoplastic 
analysis  with  the  same  defect,  spacing,  and  stress 
intensity.  We  see  in  Figure  4.30  that  the  stresses  are  in 
the  range  of  180  ksi ;  and  the  strains  are  about  six  percent 
at  the  crack  and  less  than  one  percent  at  the  defect. 

Figures  4.37-4.38  show  the  increasing  strain  profile  beivtaMin 
t.te  crack  and  defect  as  the  spacing  is  decreased.  We  note 
here  that  the  stress/'s train  data  is  comparable  to  the 
larger  primary  crack  indicating  that  the  crack  length  still 
has  little  effect  on  the  influence  range  of  the  defect.  It 
is  important  to  note  that  the  crack  stress  intensity  was 
held  constant  for  both  cases  which  required  an  increase  in 
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STRAIN  EY  (PERCENT) 

FIGURE  4.38  VISCOPLASTIC  STRESS-STRAIN  RESPONSE  DURING  CYCLIC  LOADING 


200.00 


load  for  the  second  analysis  from  18  kips  to  25  kips.  As 
long  as  similitude  is  maintained  with  the  stress  intensity* 
it  appears  the  the  defect  size  and  its  location  are  the 
primary  factors  affecting  the  material  behavior  between  the 
crack  and  defect. 

The  analysis  confirms  that  a  defect  has  a  prominent 
influence  range  equal  to  about  one  defect  length  for  all 
three  constitutive  models.  In  addition*  a  critical  spacing 
exists  in  the  plastically  deformed  region  where  the 
probability  of  instantaneous  fracture  to  the  defect  is 
imminent.  Although  no  definite  {parameter  exists  to  confirm 
the  critical  spacing*  the  critical  strain  and  growth  rate  Ca 
stress  criterion^  appear  to  be  two  primary  indicators  for 
estimating  when  instantaneous  failure  will  likely  occur. 
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A  two-dimensional  nonlinear  finite  element  code  called 
SNAP  was  used  to  model  a  center  cracked  axial  tension 
specimen  with  two  defects  placed  symmetrically  about  the 
crack  plane.  The  crack^defect  configuration  was  modeled  via 
linear  elastic,  elastic-plastic,  and  viscoplastic 
constitutive  equations.  The  material  modeled  was  IN-718,  a 
nickel  based  superalloy  used  in  Jet  engines. 

The  following  statements  and  conclusions  are  based  on 
the  analysis  presented  in  this  work: 

1.  Linear  elastic  modeling  through  finite  element 
techniques  provided  comparative  results  to  known  infinite 
plate  theory  for  a  collinear  crack/defect  configuration. 
Plots  of  the  crack  stress  intensity  versus  crack/'defect 
spacing  predict  a  critical  spacing  at  which  the  crack  stress 
intensity  will  exceed  the  combined  crack  and  defect  stress 
intensity.  This  phenomenon  has  been  referred  to  as  “pop-in‘* 
behavior  and  reveals  that  a  defect  can  actually  impede  crack 
growth  when  it  combines  with  the  primary  crack.  In 
addition,  the  defect  can  shield  the  crack  tip  from  finite 
width  effects. 

2.  A  defect  has  a  prominent  influence  range  equal  to 
approximately  one  defect  length  for  the  three  constitutive 
models  used  In  this  analysis.  Inside  this  range,  the 
defect  influence  on  the  crack  opening  displacement  and 


str«ss/^strain  flolds  cannot  be  ignored. 

3.  The  crack  sizes  considered  in  this  analysis  had  little 
affect  on  the  crack/defect  interaction.  The  primary  driver 
app>ears  to  be  the  size  of  defect  and  the  spacing  from  the 
crack  tip. 

4.  Consideration  of  the  strain  profiles  in  the  elastic- 
plastic  analysis  confirms  the  critical  spacing  concept 
revealed  with  the  linear  analysis.  Selection  of  an  eight 
percent  strain  criterion  for  failure  increases  the  critical 
spacing  over  linear  analysis  by  as  much  as  50  percent.  As 
discussed  earlier,  the  selection  of  the  critical  strain  is 
based  on  the  work  of  Mercer  CSI).  In  his  analysis,  eight 
percent  strain  was  needed  in  the  crack  tip  element  to  match 
experimental  crack  openings  for  a  compact  tension  specimen 
of  1N~718  subjected  to  a  35  KslVin  stress  intensity. 

5.  Crack  tip  blunting  and  reduced  strain  energy  at  the 
crack  tip  were  two  additional  effects  present  with  the 
nonlinear  analysis.  The  stress  intensity  calculated  within 
the  plastically  deformed  crack/defect  region  decreased  by 
approximately  50  percent  and  the  degree  of  crack  tip 
blunting  appeared  to  be  independent  of  the  crack/'defect 
spacing. 

6.  Crack  growth  with  the  viscoplastic  model  showed  that 
the  defect  will  restrict  the  plastic  zone  from  propagating 
with  the  crack  tip.  This  created  a  smaller  “super-^s trained" 
region  inside  the  original  strain  hardened  region  as  the 
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crack  continued  to  grow  towards  the  defect. 

7.  Cyclic  crack  growth  with  the  viscoplastic  model 
allowed  the  growth  rate  Cda/dtD  to  be  used  as  a  failure 
criterion.  This  method  also  predicted  a  30  percent  larger 
critical  spacing  over  the  linear  stress  intensity 

predi  ction. 

8.  The  critical  strain  and  growth  rate  criterion  are  not 
exact  methods  for  determining  failure  between  a  crack  and 
defect.  However,  they  do  provide  a  plausible  range  for  the 
problem. 

The  results  presented  in  this  work  confirm  the  ability 
of  finite  element  analysis  as  a  tool  to  imder stand 
crack^defect  interaction.  Furthermore,  some  insight  into 
how  defects  can  alter  crack  growth  parameters  under 
nonlinear  considerations  has  been  provided. 


07 


Bibliography 


1.  Military  Standard:  Engine  Structural  Integrity  Program 
(ENSIP) .  MIL-STD-1783  (USAF),  30  November,  1984. 

2.  Haritos  G. ,  Program  Manager  Aerospace  Sciences, 
AFOSR/NA,  “Structural  Durability"  Address  to  Air  Force  In¬ 
stitute  of  Technology,  Wright-Patterson  AFB  OH,  5  Nov  1986. 

3.  Zahoor  A.,  and  Abou-Sayed  I.S.,  "Prediction  of  stable 
Crack  Growth  in  Type  304  Stainless  Steel",  Computers  and 
Structures.  Vol.l3,  137-144,  1981. 

4.  Hinnerichs,  T.D.,  "Viscoplastic  and  Creep  Crack  Growth 
Analysis  by  the  Finite  Element  Method",  PhD  Dissertation, 
AFIT/DS/AA/80-2 ,  Air  Force  Institute  of  Technology,  Depart¬ 
ment  of  Aeronautics  and  Astronautics,  September,  1980. 

5.  Nicholas  T. ,  Palazotto  A.,  and  Bednarz  E. ,  "An  Analyt¬ 
ical  Investigation  of  Plasticity  Induced  Closure  Involving 
Short  Cracks",  Mechanics  of  Fatigue  Crack  Closure.  ASTM  STP 
982. pp  361-379,  1988. 

6.  Wilson  R.  and  Palazotto  A.,  "Viscoplastic  Fatigure  in  a 
Superalloy  at  elevated  Temperatures",  Fracture  Mechanics. 
Vol.  17,  ASTM  STP  905,  pp  265-275,  1986. 

7.  Palzotto  A.,  and  Henkel  C.,"The  Effects  of  Low  Cycle 
Fatigue  Comparing  Compact  Tension  and  Center-Cracked  Speci¬ 
mens  at  Elevated  Temperatures",  Engineering  Fracture  Mecha¬ 
nics.  Vol. 24,  483-494,  1986. 

8.  Mercer  d.G. ,  "Viscoplastic  Analysis  of  Fatigure  Cracks 
at  Notches  by  the  Finite  Element  Method",  PhD  Dissertation, 
AFlT/DS/AA/86-2,  Air  Force  Institute  of  Technology,  Depart¬ 
ment  of  Aeronautics  and  Astronautics,  September,  1986. 

9.  Chestnut  G.L.,  “Cyclic  Crack  Growth  Emanating  at  a 
Round  Notch  Considering  Viscoplasticity",  Master's  Thesis, 
AP1T/GAE/AA/67D-3,  Air  Force  Institue  of  Technology,  Depart 
ment  of  Aeronautics  and  Astronautics,  December,  1967. 

10.  Rice  J.R.,  "A  Path  Independent  Integral  and  the  Appro¬ 
ximate  Analysis  of  Strain  Concentration  by  Notches  and 
Cracks,"  Journal  of  Applied  Mechanics.  Vol.  35,  379-386, 
1968. 


11.  Owen  D.R.,  and  Hinton  B. ,  "Finite  Blenents  in  Plasti¬ 
city".  Pineridge  Press,  Swansea,  U.K. ,  1980. 


Bibliography  (Cont*d) 


12.  Bodner  S.R.,  and  Partom  Y.,  "Constitutive  Equations  for 

Elastic-Viscoplastic  Strain-Hardening  Materials,"  Journal  of 
Applied  Mechanics.  Vol.  42:  385-389,  1975. 

13.  Brockman,  Robert  A.,  SMAPt  A  Simple  Mon linear  Analysis 
Program  for  Education  and  Research.  University  of  Dayton 
Research  Institute  Report:  UDR-TM-82-06,  February,  1982. 

14.  Newman  J.C.,  "Finite  Element  Analysis  of  Fatigure  Crack 
Propagation  Including  the  Effects  of  Crack  Closure",  PhD 
Dissertation,  Dept,  of  Engineering  Mechanics,  Virginia  Poly¬ 
technic  Inst.,  1974. 

15.  Matake  T. ,  and  Imai  Y. ,  "Pop-In  Behavior  Induced  by 
Interaction  of  Cracks",  Engineering  Fracture  Mechanics.  Vol. 
9,  17-24,  1979. 

16.  Rose  L.R.F.,  "Microcrack  Interaction  With  a  Main 
Crack",  International  Journal  of  Fracture.  Vol.  31,  233-242, 
1986. 


17.  Chang  R. ,  "On  Crack-Crack  Interaction  and  Coalescence 
in  Fatigue". Engineering  Fracture  Mechanics.  Vol.  16,  683- 
693,  1982. 

16.  Yoda  M. •  "Subcritical  Crack  Growth  and  Crack  Coales¬ 
cence  in  Glass",  Engineering  Fracture  Mechanics.  Vol. 29, 
189-196,  1988. 

19.  Ang  W.T.,  "A  Boundary  Integral  Solution  for  the  Problem 
of  Multiple  Interacting  Cracks  in  an  Elastic  Material", 
International  Journal  of  Fracture.  Vol.  31,  259-270,  1986. 

20.  Rubinstein  A. A,,  "Hacrocrack  -  Hicrodefect  Interac¬ 
tion",  Journal  of  Applied  Mechanics.  Vol.  53,  505-510,  1986. 

21.  Yanada  Y.,  et  al.,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  for  the  Solution  of  Elastic-Plastic  Problems 
by  the  Finite  Element  Method". International  Journal  for 
Mechanical  Structures.  Vol.  10,  343-354,  1968. 

22.  Hughes  T.J.,  and  Taylor,  "Unconditionally  Stable  Algo¬ 
rithms  for  Quasi-Static  Blasto/Viscoplatic  Finite  Element 
Anal  vs is ".Computer a  and  Structures. Vol .  0,  Ho. 2,  169-173, 
1978. 


23.  Beaman  R.,  "The  Determination  of  the  Bodner  Material 
Coefficients  for  XH-716  and  Their  Effects  on  Cyclic  Load¬ 
ing",  M.S.  Thesis,  AF1T/GAB/AA/64M-1 ,  School  of  Engineering, 
Air  Force  Institute  of  Technology,  Wright-Patterson  AFB, 
Ohio,  1983. 


Bibliography  (Cont*d) 


24.  Broek,  D.,  Elementary  Fracture  Mechanics.  (Fourth  Edi¬ 
tion),  The  Hague:  Martinus  Nijhoff  Publishers,  1986. 

25.  Isida  M. ,  "Effect  of  Width  and  Length  on  Stress  Inten¬ 
sity  Factors  Of  Internally  Cracked  Plates  Under  Various 
Boundary  Conditions",  International  Journal  of  Fracture 
Mechanics.  Vol.7,  301-316,  1971. 


26.  Cook  R.D.,  Concepts  and  Applications  of  Finite 
Element  Analysis.  (Second  Edition),  John  Wiley  &  Sons,  1981. 


Appendix  A 

Simple  Nonlinear  Analysis  Program 
SNAP 

The  following  background  and  user  information  is 
provided  to  reflect  the  changes  made  by  Mercer  C8^  and  Smith 
to  the  AFIT  version  of  the  SNAP  Finite  Element  Code  Cl 35. 

The  source  code  title  is  SNAP16D.  F  and  is  compatible  with 
all  VAX  11-75  compilers. 

GENERAL  DESCRIPTION 

SNAP  is  a  finite  element  solution  program  for 
two-dimensional  problems  in  nonlinear  structural  mechanics. 
The  program  can  be  used  foi'  problems  with  the  following 
character i St i cs : 

-  plane  stress  ot  plane  strain  assumptions 
~  small  or  moderately  large  displacements 

-  arbitrary,  two-dimensional  geometry 

-  cyclic  load  y  monotonic  loading 

-  crack  propagation*  including  opening  and  closure 

-  model  size  less  than  6000  O.  O.  F. 

-  elastic  or  time-independent  elastic -plastic  materials 

-  viscoplastic  or  time-dependent  analysis 


DATA  PROCESSING  ft  FILE  gTRUCTURE 

Data  processing  Is  kept  in  core  under  control  of  a 


simple  dynamic  storage  allocation  scheme;  the  only  use  of 
secondary  storage  is  in  the  maintenance  of  element  state 
variables  and  the  creation  output  files  requested  by  the 
user. 

Operation  of  the  snap  program  requires  the  following 
files: 


File 


Description 


snapin.  dat 
snapout . dat 
di spl , dat 
sigstra. dat 
r St new 
rstold 
time 


primary  data  input  deck 
primary  data  output  deck 
node  displacement  output  file 
element  stress/strain  output  file 
restart  file  created  at  end  of  run 
existing  restart  file  from  prior  run 
dump  file  for  processing  status 


INPUT  DATA  C€^VENTtO»<S 

Problem  input  for  SNAP  is  divided  into  several  **data 
blocks",  which  are  used  as  needed  to  define  the  problem  to 
be  solved.  A  summary  of  the  available  data  blocks,  their 
functions,  and  whether  or  not  they  are  required  for  all 


problems  is  given  below. 

Block  Description  Requi  red 


B01A4  Define  nodal  constraints  Y£S 

COOR  Define  nodal  coordinates  YDS 

CRAG  ^Mtcify  nodes  and  directions  for  ND 

crack  opening  simulation 


BLOCK  INSCRIPTION  REQUIRED 

CLOS  Sptecify  nodes  and  directions  for  crack  NO 

closure  simulation 

ELEM  Define  finite  elements  YES 

FORC  specify  concentrated  forces  at  nodes  NO 

HIST  Define  loading  as  a  function  of  time  YES 

JINT  Define  J-integral  paths  NO 

MESH  Generate  a  scaled  plot  of  elements  NO 

PARA  Select  type  of  analysis  and  solution  YES 

parameters 

PRES  Specify  distributed  element  loads  NO 

PRIN  Print  specific  increments  of  load  history  NO 

PROP  Define  elastic  material  properties  YES 

REST  Create  a  restart  file  for  a  problem  NO 

STAT  Initialise  material  state  variables  YES 

SOLN  Solution  start /stop  after  input  complete  NO 

SFIL  Define  elements  for  stress  and  displace^  NO 

ment  files  Csigstra.dat  &  displ.dat > 

TITL  Define  descriptive  title  for  a  problem  YES 

INPUT  FORMAT 

Each  input  data  block  begins  with  a  single  line 
containing  the  header  of  the  block  Ce«  g.  '^OOOR’O.  Data 
blocks  may  be  entered  in  any  CM*der»  but  each  block  may 
appear  only  once  in  the  input  file. 

BOUNDARY  CONDITIONS  Input  Block 


Headei^:  BOUN 


BOUNDARY  CC»4D1TI0NS  Input  Cont’d 
End:  Blank  line 

Record  Format  C repeat  as  needed^ 

. 5. _ 10 . 15 . 20 . 25 

IBEG  lEND  INCR  IXFIX  lYFIX 

Variables: 

IBEG  =  Beginning  node  to  be  constrained 
BOUNDARY  CONDITIONS  Input  Cont»d 

I END  ~  Ending  node  to  be  constrained 
INCR  =  Node  number  increment 

IXFIX  =  X'direction  constraint  code  CO  =  free»  1  =  fixed^ 
lYFIX  »  Y -direct! on  constraint  code  CO  =  free,  1  »  flxedD 

COORDINATES  Input  Block 

Header:  COOR 

END:  blank  line 

Record  Format  CID: 

...  .5 
NUMNOD 

Variables: 

N18<INCX>  s  Total  number  of  nodes  def i  ned  for  the  ncxlel . 

Record  Format  C23: 

.  .5 . 10 . 20 . 30 

NODE  I  GEN  XCORD  YCORD 

Variables: 

NODE  »  Node  number 

I  GEN  s  Generation  increiment 

XCORD  «  X-coordinat« 


YCORD  *  Y -coordinate 


CRACK  Input  Block 


Hoader :  CRAC 

End:  CnoheD 

Racord  Format  Cl^: 

. B 

NUMCON 

Variable: 

NUMCON  s  Number  of  conditional  nodes  in  crack  front 

Record  Format  C8^: 

_ S . 10 

NODE  XDIR 
Variables: 

NODE  3  Node  number  for  crack --opening  constraint 

IDXR  «  e  Direction  2«V5  of  constraint 

CM  positive  value  closed*  negative  openD 

Record  Format  C33: 

. 6 . 10 . IS 

KRKTfPE  NSTTEP  IPOP 

Variables: 

KRKTYPE  a  Type  of  element  used  In  crack  front 

2  **  Four  noded  quadr  lateral 

3  -  Eight  noded  quadrilateral 

NSTEP  «  Number  of  increments  to  reach  full  load 

I POP  s  Number  of  nodes  released  per  cycle 

Record  Format  CO: 

. 6 

1  STRESS 


CRACK  Input  Cont'd 
Variable: 

ISTRESS  s  Stress  required  to  release  a  node 

CRACK  CLOSURE  Input  Block 

Header :  CLQS 

End:  CnoneD 

Record  Format  Cl^: 

. 5 

NUMCON 

Variable: 

NUMCON  «  Number  of  conditional  nodes  in  crack  front 

Record  Format  CS>t 

. S 

SOPEN 

Variable: 

SOPEN  »  Stress  to  c^n  closed  nodes 
Record  Format  C33: 

- S . 10 

NC£)£  IDIR 
Variables: 

NODE  s  Node  number  for  crack '-opening  constraint 

IDXR  s  M  Direction  Ci=X.  of  constraint 

Cm  initially  all  nodes  have  a  negative  valued 

ELEMENTS  Input  Block 
Header:  ELEN 


.End:  blank  line 


ELEMENTS  Input  Cont’d 
Record  Format  Cl  I): 

...  .5 
NELEM 
Variable: 

NELEM  -  Number  of  elements  in  model 
Record  Format  C2D: 

. 6 . 10 . 15 . 30 . 25 . 30 . 36 

lELEM  lEGEN  INGEN  NODECID  NODEC23  N0DEC3i  NODEC45 

. 40 . 45 . 50 . 55 . 60 

N0DEC53  N0DEC6D  NODEC75  N0DEC85  INTORD 

Varlabl  es: 

lELEM  s  Element  number 

I GEN  s  Element  generation  increment  C»  not  usedD 

INGEN  =  Node  number  increment  generation  C**  not  used^ 

NODE  CiD  s  i~th  connected  node  for  this  element 

INTCKD  =  Numerical  integration  order 

Cl  or  3  =  1 -point  or  3-point  rule,  for  triangles^ 
C4  or  Q  s  2x2  or  3x3  Gauss  rule,  for  quadrilat¬ 
eral  sD 

M  lEGEN  AND  INC^N  are  not  used  if  mesh  is  generated  with 
external  program  such  as  PDA-PATRAN. 


Forces  Input  Block 
Header:  FORC 

End:  blank  line 

Record  Format:  C repeat  as  needed!) 


.  .  .5 . 10 . 20 

NODE  I DIR  FORCE 


Variables: 
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FORCES  Input  Cont'd 

NODE  -  Node  number  at  which  the  force  is  applied. 

IDIR  =  Direction  Cl  =  X,  2  =  ¥5  of  the  force 
FC®CE  =  Force  magnitude 

HISTORY  Input  Block 

Header :  HI  ST 

End;  CnoneD 

Record  Format  Cl!): 

. 5 

NUMPTS 

Variable: 

NUMPTS  =  Number  of  p>oints  Ct.  fCt^^  to  be  defined  on  force 
time  history  curve 

Record  Format  CSD:  C repeat  as  needed^ 

. 10 . 20 

TIME  FUNG 

Variabl es: 

TIME  s  Reference  value  of  time  Cloadlng  parameter^ 

FUNG  B  Corresponding  scale  factor  for  applied  loads 

J -INTEGRAL  Input  Block 

Header:  JINT 

End;  Cnone^ 

Record  Format  CID: 

. 20 

TITLE 

Variable: 

TITLE  •  20  character  title  block 


J -INTEGRAL  Input  Cont*d 


Record  Format  C23: 

. 5 . 10 

INUMBER  MCX)EL 

Variables: 

INUMBER  =  Number  of  Paths 

MODEL  =  Type  of  Constitutive  Model 

1  —  Elastic 

2  —  Elastic  Plastic 

Record  Format  C3D: 

...  .5 
ETOTAL 

Variable: 

ETOTAL  =  Total  Number  of  Elements  in  path 

Record  Format  C43:  C repeat  as  needed!) 

. 5 . 10 . IS . 20 . 2S . 80 

ELEMC 10  ELEMC  20  ELEMC  30  ELEMC  40  ELEMC  5^ .  .  .  ELEMC 1 6D 

Variable: 

ELEMCiD  s  i-th  element  number  containing  J-integral  path 

Record  Format  C53:  Crepeat  as  neededO 

. 8 . 10 . IB . 20 . as . 80 

IDIRC13  1DIRC20  IDIRCa:)  IDIRC40  IDIRCB^ .  .  .  IDIRCIOO 

Variable: 

IDIRCi!)  a  1-th  path  direction  CSee  Diagram  BelowD 
e  NOTE  EACH  ELEMENT  HAS  TWO  DIRECTIONS  *• 


ELEMENT  PATH  DIAGRAM 


•  dimkcticmn  ■  ur  mhd  to  tmk  micht 

->  eiKKCTIOM  m  eOMM  NMO  YO  VMIC  UKVY 
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MESH  PLOTTING  Input  Block 
Header:  MESH 
End:  Cnone^ 


Record  Format  CID: 

. 5 . 10.  .....  20.  ....  30.  ....  40 . 50.  .....  CO 

lOPT  INUM  SCALE  XMIN  XMAX  YMIN  YMAX 

Variables: 

I^T  s  Option  for  mesh  plotting 
O  —  No  plot 
1  —  Plot  mesh 

INUM  =  Option  for  plotting  element  numbers 
O  —  No  numbers 
1  —  Plot  numbers 

SCALE  3  Plotting  magnification  scale 
1-OQ  —  allowable  values 

XMIN  a  Minimum  X -value  of  plot  region 

XMAX  a  Maximum  X-value  of  plot  region 

YMIN  a  Minimum  Y-value  of  plot  region 

YMAX  a  Maximum  Y-value  of  plot  region 

mm  This  routine  is  disabled  in  snaplCd.f  since  no  1038 
device  driver  is  available  on  the  AFXT/ICC.  ASD/CYBCR  ver¬ 
sion  of  SNAP  CsnaplCc.f^  supports  this  routine. 


PARAMETERS  Input  Block 
Header:  PARA 
End:  CnoneS 

Record  Format  CID: 

. . 5 . 10 . .  15. . . .  .  30.  . 

IPLANE  XSMALL  XMOOEL  NCSTIEP 

Variables: 

IPLANE  «  Problem  type 

0  —  Plane  strain 
i  —  Plane  stress 


.35. 


.  30 


XSPRT  IGPT 
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PARAMETERS  Input  Cont»d 

I  SMALL  =  Kinomatics  option 

0  —  Large  displacements 
1  —  Small  displacements 

IMCX3EL  =  Constitutive  model  type 

1  —  Elastic 

2  —  Elastic,  perfectly  plastic 
5  —  Bodner-Partom  viscoplastic 

10  —  Bilinear  Elastic  plastic 

NCSTEP  s  Total  number  of  tine  step  changes 

ISPRT  =  integration  point  stress^'strain  data  flag 
O  —  none 
1  —  Print  data 

IGPT  =  Gauss  point  coordinate  flag 
O  —  no  data  print 
1  —  print  coordinate  data 


Record  Format  C2D: 


_ 5 . 10. 

NSTEP  IPFREQ 


,  .  .  IS . 20 . 30,  ......  40 . BO 


ITOPT  MAXIT  EQTOL  WSTOL  DELTAT 


Variables: 


NSTEP  =  Number  of  solution  increments 


IPFREQ  =  Printing/file  output  frequency  Cin  increments!) 

ITOPT  a  Nonlinear  solution  option 

0  —  Corrected  step-by-step 

1  —  Constant  stiffness  C pseudof orce> 

2  --  Modified  Newton  iteration 

3  —  Full  Newton  Iteration 

4  —  Combined  Newton  Iterations 

MAXIT  a  Maximum  number  of  iterations  per  step 

EQTOL  a  Convergence  tolerance  on  force  residuals 

DISTOL  a  Convergence  tolerance  on  displacement  corrections 

DELTAT  «  Time  Cload  parameter?  step  size 
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PARAMETERS  Input  Cont*d 

Record  Format  CSD  C repeat  as  needed^ 

. 5 . 15 .  20 .  30 .  35 .  45 _ 80 

ICSTEPCli  CDTCID  ICSTEPC23  CDTC2i  ICSTEPCS^  CDTC3D . 

CMAXIMUM  STEP  CHANGES  -  2003 

Variables: 

ICSTEP  Ci3  s  1— th  increment  or  step 
CDTCi3  =  i-th  time  step  size 

This  record  enables  the  user  to  apply  different  time  step 
increments  during  a  load  cycle. 

PRESSURES  Input  Block 

Header:  PRES 

End:  C blank  lineD 

Record  Format:  Crepeat  as  needed^ 


.  .5 . 10 .  80 .  30 .  40 

lEL  lEOGB  PRESSA  PRESS8  PRESSC 


Variables: 

lEL  «  Element  to  which  pressure  is  applied 
lEOGE  «  Edge  Number 

1  B  top.  2  B  left  side*  3  «  bottom  4  •  right  side 
PRESSA  B  Pressure  value  at  first  node  point  on  the  edge 
PRESS8  B  Pressure  value  at  second  node  point  on  the  edge 
PRESSC  B  Pressure  value  at  third  node  point  on  the  edge 

PRINT  Input  Block 
Header:  PRIN 
End:  Cnone) 
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PRINT  Input  ContM 
Record  Format  Cl^: 

. 5 

I  PRINT 

Variable: 

IPRINT  =  Number  of  increments  to  print  CMAXIMUM  ROOD 

Allows  user  to  print  increment  that  are  skipped  by 
IPFREO  in  PARAMETERS  block. 

Record  Format  C2^: 

. 5 . 10 . 15 .  20 .  25 .  80 

ICSTEPCID  ICSTEPC2D . 

Variable: 

ICSTEPCi^  =  i-th  increment 

PROPERTIES  Input  Block 

Header:  PROP 

End;  Cnone!> 

Record  Format: 

. 10 . 20 

ELAS  POIS 

Variables: 

ELAS  s  Elastic  modulus 

POIS  «  Poison's  ratio  0  <  POIS  <0.50 

STATE  VARIABLES  Input  Block 
Header:  STAT 

End:  CnoneD 

Record  Format  CIO: 

. 6 . 10 

NSTATE  NCONST 


Variables: 


STATE  VARIABLES  Input  Cont*d 

NSTATE  =  Number  of  Internal  state  variables  to  be  defined 
NCONST  s  Number  of  constant  parameters  to  be  defined 

Record  Format  C2D:  C repeat  as  needed^ 

. 10 . 20 . 30 . 40 . 80 

STATECIO  STATEC2D  STATE  C 3D  STATECAD.  .  .  STATE  C8D 

Variables; 

STATECiD  =  Initial  value  of  the  i-th  internal  state 
variable. 

Record  Format  C3D:  Crepeat  as  neededD 

. 10 .  20 .  30 .  40 .  80 

CONSTCID  C0NSTC2D  CONST  C 3D  CONSTCAD.  .  .  CONST  C8D 

Variable: 

CONSTCiD  s  Value  of  i-th  constant  parameter  value. 

C  MAXI  MUM  20D 

STTRESSXSTRAIN  &  DISPLACEMENT  FILES  Input  Block 
Header :  SFI L 
End:  CnoneD 

Record  Format  CID:  Crepeat  as  neededD 

. S . 10 

lELEMClD  ISIRCID 

Variables: 

lELEM  Number  of  elements  to  be  output  C200  maximumD 

ISTR  a  Stress/'strain  Designator 

1  »  XX 

2  »  YY  C DEFAULT) 

Record  Format  C2D:  Crepeat  as  neededD 

. 10 .  20 .  30 .  40 .  80 

JELMCID  JELMC2D  JELM  C3D  JELNC4D . JELNC8D 

Variables: 
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STRESS/STRAIN  &  DISPLACEMENT  FILES  Input  Cont'd 


JELMC13  =  Eloment  numbers  whose  stress/strain  will  be  written 
sigstra.  dat 


Record  Format  C3^: 

. 5 

NNODES 

Variables: 

MNOIi^IS  =  Number  of  Nodes  whose  displacement  will  be  written 
to  di  spl .  dat 


Record  Format  C43: 

. 6 . 10 . 15 .  20 .  25 .  80 

MNODECID  IDIRCID . 

Variabl es: 

MNODECi^  -  i-th  node  number 

IDlRCiD  8  i>th  displacement  direction 
1  8  X 
a  8  Y 


TITLE  Input  Block 
Header:  TITL 

End:  CnoneD 

Record  Format: 


80 


Variables: 

TITLE  8  Alphanumeric  problem  title  Cup  to  80  characters^ 


MATERIAL  MODEL  VARIABLE  DESCR1PTICX4S 

XMOOEL  -  1:  LINEAR  ELASTIC  MATERIAL 
C^stants:  Cnone!) 


State  Variables:  CnoneS 


I  MODEL  =  2;  ELASHC  -  PERFECTLY  PLASTIC 
Constants:  CnoneD 

State  Variables: 

1.  Yield  Stress 

2.  Effective  Stress  Clnitially  set  to  zero^ 

3.  Elastio^lastic  Plastic  Work 


I MODEL  =  S: 

BODNER-PARTOM  VISCOPLASTTIC  MATERIAL 

Constants: 

1. 

DO 

2. 

ZO 

3. 

21 

4. 

22 

S. 

n  Cfloating  point!? 

6. 

m  C floating  point? 

7. 

A  Cfloating  point? 

State  Variabl 

es: 

1. 

Plastic  Strain  Exx 

2. 

Plastic  strain  Eyy 

3. 

Plastic  strain  Ezz 

4. 

Plastic  strain  2*^b(y 

B. 

2 

0. 

Plastic  work 

7. 

Total  Work  CElastlc  ■*>  Plastic? 

INOOEL  e  10:  EIUNEAR  ELASnC-PLASTlC  MATERIAL 


104 


Constants: 


1.  Yiold  strass 


2.  Strain-’hardaning  slope 
State  Variables: 

1.  Plastic  Strain  In  X-dlrectlon 

2.  Plastic  Strain  In  Y -direction 

3.  Plastic  Strain  in  Z-dlrectlon 

4.  Plastic  Strain  in  XY-dlrectlon 

6.  Effective  Plastic  Strain 

6.  Effective  Stress 

7.  Current  Yield  Stress 

8.  Elastic  Plastic  Work 
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APPENDI X  B 


SWAP  PROGRAM  SOURCE  CODE 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 
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PROGRAM  SWAP  C INPUT, OUTPUT, CHOICES  VIA  SCREENS 

THIS  PROGRAM  ACTS  AS  AN  INTERFACE  BETWEEN  THE  PATRAN 
NEUTRAL  FILE,  THE  NASTRAN  BANDIT  ALGORITHM,  AND  THE 
SNAP  INPUT  DATA  FILE 

THREE  OPTIONS  ARE  AVAILABLE  TO  THE  USER: 


1.  GENERATE  BANDIT  INPUT  DATA  FILE  FROM 
THE  PATRAN  NEUTRAL  FILE 

a.  GENERATE  A  SNAP  INPUT  DECK  USING  THE 
USING  THE  PATRAN  NEUTRAL  FILE  AND 
AND  BANDIT  RESULTS 

3.  GENERATE  SNAP  INPUT  DATA  FROM 
PATRAN  NEUTRAL  FILE  ONLY 


FILES  ARE  AS  FOLLOWS: 

PATRAN.DAT  ~  PATRAN  NEUTRAL  FILE  CINPUD 
BANDIT.DAT  —  BANDIT  INPUT  DATA  C NASTRAN  OUTPUD 
BANDIT.  LOG  —  BANDIT  RENUMBERING  RESULTS  CINPUD 
SNAP.INP  --  SNAP  INPUT  DATA  COUTPUD 

DI MENSI ON  NCONC  3000 , 80 , XC  40000 , YC  40000 
DI MENSI ON  I OLDC  40000 , 1 NEWC 40000 
CHARACTER'S  I TYPE 

OPEN  CUNIT«8.FILE='PATRAN.DAT» ,STATUS=»OLD’0 
OPEN  C  UNI  T==0 ,  FI  LE= » BANDI T.  DAT’  ,  SrATUS= » NEW*  0 
OPEN  C  UNI T=1 0 , FI LE=  *  BANDI T.  LOG’ , STATUS= ’ OLD*  0 
OPEN  C UNI Tal 1 , FILE* ’SNAP.  INP’ , STATUS* » NEW ’0 

REWIND  8 
REWIND  0 
REWIND  10 
REWIND  11 
2  *  0.0 
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n  PROGRAM  SWAP  CONT'D  ** 

WRITEC6»aOOOD 
WRITEC6,ai00D 
READ  C5,»tD  lOPT 

IF  CIOPT  .EQ.  1  ?  WRITE  C9,5000D 
IF  CCIOPT  .NE.  2D  .AND.  CIOPT.  NE.  3DD  GO  TO  50 
WRITEC6,1050D 
READC5,»*D  INTORD 
50  CC»ITINUE 

C 

READ  C 8,10005  NOCCS.NEL 
CC  WRITEC6.*tD  NOi:*:S,NEL 

DO  1 00  I  =  1 , NODES 
READ  C 8, 9605  INODE 
READ  C 8, 11005  XCI N0DE5 , YCI N0DE5 
READ  C 8, 9605 

CC  WRITE  C6,it5  INODE 

CC  WRITE  C6,*<5  XCI NCOED  ,  YCI NODE5 

100  WRITE  C9, 11505  IN0DE.XCIN0DE5  ,  YCIN0DE5 ,2 

C 

C  READ  ELEMENT  CONNECTIVITY  Fe*OM  PATRAN  NEUTRAL  FILE 

C 

DO  200  I  =  1,NEL 
READ  C 8. 9705  I EL 
READC8, 12005  CNCONCIEL,  J5  ,  J-1 ,85 
CC  WRITE  C6,»«5  lEL 

CC  WRITE  C6,»«5  CNCONCIEL,  J5  ,  J*1 ,85 

ITYPE  =  ’cisaoe  • 

ITYP2  =  *+CIS  • 

MID  =  1 
ITEST  =  O 

IF  CNC0NCIEL,55  .  GT.  05  GO  TO  22 
ITfPE  =  ’CGUADl  ’ 

ITEST  =  1 

22  IF  CNC0NCIEL.45  .  GT.  05  GO  TO  23 

ITYPE  =  'CTRI 

ITEST  =  2 

23  IF  CIOPT  .GT.  15  GO  TO  200 
IF  CITEST  .EQ.  15  GO  TO  24 
IF  CITEST  .EQ.  25  GO  TO  25 

WRITE  CO, 22005  ITYPE,IEL,MID,CNCONCIEL,  J5  ,  J=1 ,05  , 
ITYP2,IEL,  ITYP2,IEL,CNC0NC1EL,K5,K«7,85 
GO  TO  200 

24  WRITE  CO,  22245  ITYPE, lEL, MID,CNCONCIEL,  J5  ,J*l  ,45 
GO  TO  200 

25  WRITE  CO, 22255  ITYPE, lEL, M1D,CNC0NCIEL,  J5  ,J»1 ,35 
200  CONTINUE 

C  END  IP  GENERATING  NASTRAN  TYPE  DATA  FOR  BANDIT  INPUT 
C 

IF  CIOPT  .EQ.  15  WRITE  CO, 32005 
IF  CIC3PT  .EQ.  15  WRITECO, 60105 
IF  CIOPT  .EQ.  15  GO  TO  000 
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«  PROGRAM  SWAP  CONT*D 

IF  CIOPT  .EQ.  3D  GO  TO  710 

C  READ  RENUMBERING  DATA  FROM  BANDIT  OUTPUT 
C 

READ  C10,1500D  CIOLDC  JD  »INEWC  JD  ,  J=1  .NODESD 
C  DO  300  I =1 , NCOES 

C  WRITEC6,*»D  lOLDCID.INEWCID 

C  WRITEC6,**D  XCID»YCID 

C300  C(»ITINUE 
C 

C  CONVERT  NODAL  COORDINATES 

C 

DO  400  I  =  1, NODES 
M  =  I 

N  =  INEWCID 
XCND  =  XCMD 
YCND  =  YCMD 

C  WRITEC6,»*D  M.INEWCID 

400  CONTI  HUE 
C 

C  CONVERT  ELEMENT  CONNECTIVITIES 

C 

DO  SOO  lEL  =  l.NEL 
C 

DO  450  J  =  1 .8 
M  «  NCONCIEL.JD 
IF  CM  .  LE.  OD  GO  TO  460 
NCONCIEL.JD  a  INEWCMD 
450  CONTINUE 
C 

500  CCJNTINUE 
C 

C  WRITE  OUT  RESULTS  IN  SNAP  FORMAT 
C 

710  WRITE  Cll  ,37005 

WRITE  Cll, 27505  NODES 
DO  760  INODE  =  1, NODES 

750  WRITE  Cll, 27805  1N0DE,XCXN0DE5 , YCINODED 
WRI TECH  ,20005 
WRZ TEC 11, 28005 
WRI TEC 11, 27505  NEL 
DO  800  lEL  a  1,NEL 

800  WRI TEC 11, 28505  IEL,CNCONCIBL. J5 , J»1 ,85  ,  INTORD 
WRITE  Cll, 20005 
WRI  TEC  5,31005 
000  CONTINUE 
STOP 

060  PORMATCSX,I65 

070  FORMATCBX,I6,/5 

1000  F0RMATC^,/,30X,X5,2X,IB./5 

1060  F0RMATC41H  INPUT  ELEMENT  INTEGRATION  ORDER  FOR  SNAP  5 
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**  PROGRAM  SWAP  COIT’D  *» 

1 1 00  FORMATC  3E1 6.  OD 

1150  FORMATC 8HGRID  ,I8,8X,3F8.  5!) 

1200  FORMATC 8C4X, I 4DD 
15CX)  FORMATC  32X.  81 83 

20CX>  FORMATC  ✓/AOHXiOOtWIOtWWWWKWWMWMIOtMMMWiWKIOOtWMWMKMMiWICMX 


2060 


2100 


2200 

2224 

2225 
2700 
2750 
2780 
2800 
2850 
2000 
3100 


3200 

C 

C  «• 

C 

5000 


6010 


+ 

40H 

NODAL  RENUMBERING  PROGRAM 

FCSMATC 

40H 

^  9 1 1  uipifu  ^  0 1  j  ui 

FILE  NUMBERING: 

-f 

40H 

PATRAN.DAT  =  PATRAN  NEUT  FILE  CINP3 

4- 

40H 

BANDIT.DAT  =  DATA  FILE  BANDI TC 0UTPT3 

+ 

40H 

BANDIT.LOG  =  BANDIT  RENUMBERING  RESULTS 

+ 

40H 

CINPUT3 

+ 

40H 

SNAP.INP  =  COORDS  AND  CONNECTIVITIES 

+ 

40H 

IN  SNAP  FORMAT  C0UTPUT3 

./3 

FORMATC 

40H 

CHCX3SE  OPTIONS: 

+ 

40H 

1.  GENERATE  BANDIT  INPUT  DATA 

4 

40H 

FROM  NEUTRAL  FILE 

+ 

40H 

2.  RENUMBER  USING  BANDIT  RESULTS 

+ 

40H 

GENERATE  SNAP  DATA  FILE 

+ 

40H 

3.  REFORMAT  PATRAN  DIRECTLY  TO  SNAP 

3 

FORMATC A8 . 81 8 » A4 , I 4 . /A4 , I 4 , 21 83 

FORMATC  A8, 6183 

FORMATC  A8, 51 83 

FORMATC 4HC00R3 

FORMATC 153 

FORMATC 1 5 , 8X , aFl 0.  53 
FORMATC 4HELEM3 
FORMATC 15,1 OX, 01 53 
FORMATC 1X3 


FORMAT  C40H  SNAP  FILE  WRITTEN. 

+  40H  OUTPUT  FILE  IS:  ,</ 

+  40H  SNAP.INP  SNAP  DATA  FILE  ,/3 

FORMATC 40H  BANDIT  DATA  CBANDIT.0AT3  IS  READY  ,/3 


NOTE  CHANGE  MATERIAL  CARD  TO  NATCH  REQUIRED  ELEMENT 


FORMATC  41 HNASTRAN  TI  TLE0PT«-1 ,  BANDIPCH-l ,  BANDTRUNel ,  X 


4 

30HID  SNAP. BANDIT 

ON  MESH 

4 

30HAPP  DISPLACEMENT 

4 

30HTXNE  2 

•  ^ 

4 

30HS0L  1.0 

4 

30HCEND 

4 

30HSPC-1 

4 

30HSTRES£CS»ALL 

4 

30HBEGIN  BULK 

4 

60HSPC1 

1  345126 

1 

14 

FORMATC  60HPC»;AD1 

1  1 

.6 

4 

60HMAT1 

1  10.  ose 

.30 

4 

7HENDDATA  3 

.100 


END 


3 


114 


VITA 


LeRoy  K.  Smith  was  born  on  27  December  1960  at  Pier re »  South 
Dakota.  He  graduated  Trom  high  school  in  Mur  do »  SD  in  1070. 
He  then  attended  South  Dakota  State  University  in  Brookings » 
SD.  In  December  1083  he  was  awarded  a  B.  S.  in  Mechanical 
Engineering  with  high  honor  and  was  commissioned  as  a  second 
lieutenant  in  the  United  State  Air  Force.  He  was  then 
assigned  to  HQ  ^pace  Division.  Los  Angeles.  CA.  During  this 
assignment  he  served  as  a  facility  integration  engineer  for 
the  Automated  Remote  Tracking  Station  CARTS)  System  Program 
Office.  He  then  entered  the  School  of  Engineering.  Air 
Force  Institute  of  Technology,  in  June  1087. 


115 


UNCLASSIFIED 


EiiURITV  ClA^SIFICATION  Of  THIS  PAgf 


1  REPORT  DOCUMENTATION  PAGE  omNo%7o4%iB8 

1  U,  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb  RESTRICTIVE  MARKINGS 

None 

1  2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
Distribution  unlimited 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

1 

1  4  PERFORMING  ORGANI2ATION  REPORT  NUMBER(S) 

AFIT/GAE/AA/88D-35 

1 _ _ _ _ 

5.  MONITORING  ORGANIZATION  REPORT  NUM8ER(S) 

1  Sa  NAME  OF  PERFORMING  ORGANUATION 

School  of  Engineering 

6b.  OFFICE  SYMBOL 
(If  applicable) 

AFIT/ENY 

7a.  NAME  OF  MONITORING  ORGANIZATION 

1  ec  ADDRESS  (City,  State,  and  HP  Code) 

Air  Force  Institute  of  Technology 
Wright-Patterson  AFB  OH  45433-6583 

7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

1  ba  name  of  funding /sponsoring 

ORGANIZATION 

^\ir  Force  Wright  Aero  Lab 

8b.  OFFICE  SYMBOL 
(If  applicable) 

MLLN 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

1  8c  ADDRESS  (City,  State,  and  HP  Code) 

'  AFWAL/MLLN 

Wright-Patterson  AFB  OH  45433-6583 

1 _ 

10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO. 

PROJECT 

NO. 

TASK 

NO. 

WORK  UNIT 
ACCESSION  NO. 

■  n.  TITLE  (ln</udt  S€<uhty  CUsstfication) 

.  See  box  19 

12  PERSONAL  AUTH0R(S) 


j  LeRoy  K.  Smith 

,  Capt,  USAF 

I  I3a.  TYPE  OF  REPORT 

13b.  TIME  COVERED 

14.  DATE  OF  REPORT  (Year.  Month.  Day) 

15.  PAGE  COUNT 

1  Thesis 

FROM  TO 

1988  December 

115 

j  *  16.  SUPPLEMENTARY  NOTATION 

1 _ 

_ _ 

!T  a  O 

17. 


COSATI  COOES 


FIELD 

GROUP 

SUB-GROUP 

1 

01 

1 

SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Crack  Growth,'  Finite  Elements,'  Viscoplasticityy 

LjUcP  )  -tS= — 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

TITLE:  MACROCRACK- ^!ULTIPLE  DEFECT  INTERACTION  CONSIDERING 
ELASTIC,  PLASTIC,  AND  VISCOPLASTIC  EFFECTS 

Thesis  Chairman:  Dr  Anthony  Palazotto 

Professor  of  Aeronautics  and  Astronautics 


\v 


20  OISTRIBUnON/ AVAILABILITY  OF  ABSTRACT 
HuNClASSIFIEO/UflUMITEO  □  SAME  AS  RPT 


□  OTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22a  riAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr  A  Palazotto,  Professor 


22b.  TELEPHONE  (Include  Area  Code) 

513-255-2998 


22c.  OFFICE  SYMBOL 

AFIT/ENY 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


UNCLASSIFIED 


BOX  19: 

A  finite  element  investigation  was  conducted  to  analyze  an 
axial  tension  specimen  with  collinear  defects  placed  symmetri¬ 
cally  about  a  center  crack.  The  material  modeled  was  IN-718,  a 
nickel-based  superalloy  used  in  jet  engines.  The  effects  of 
crack/defect  interaction  were  compared  using  elastic,  elastic 
plastic,  and  viscoplastic  constitutive  models.  A  2-D  nonlinear 
finite  element  code  called  SNAP  was  used.  This  program  has  the 
capability  to  simulate  crack  growth  and  closure  by  releasing  or 
closing  nodes  along  the  crack  plane. 

Elastic  stress  intensity  solutions  were  developed  for  two 
different  finite  width  specimens.  The  stress  intensity  versus 
crack  length  plots  compared  well  with  infinite  theory.  Results 
reflect  the  defect  can  partially  shield  the  crack  from  finite 
width  effects.  A  critical  spacing  was  also  noted  where  the 
stress  intensity  of  the  crack  exceeded  the  stress  intensity  for 
the  combined  length  of  the  crack  and  defect.") 

Finite  element  analysis  of  a  crack/djjtect  configuration, 
considering  elastic-plastic  and  elaatrtS^^^scoplastic  effects, 
provided  crack  opening  prof iles, /plastic  zone  profiles,  and 
stress/strain  fields,  in  general defec.t  has  a.  prominent 
influence  range  equal  to  approximately  one  defect  length  for  all 
constitutive  models.  The  presence  of  a  defect  increases  the 
magnitude  of  the  crack  opening  and  stress/strain  fields  in  front 
of  the  crack  tip.  ^  /  /  .  . 

<^u>ov(Ur,  ^  fZ 


UNCLASS I PIE 


